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Table | 
Systems containing CO, and medium size components 











2nd Component T [K] range P [bar] range LCVM 





Ay x 1000 AP% Ay x 1000 
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Propane 274-344 
n-Heptane 310-394 
n-Octane 313-383 
n-Decane 344-377 
n-Hexadecane 463-664 
Methylcyclohexane 311-477 
Toluene 311-477 
Ethylbenzene 308-328 


Butylbenzene 273-293 
Decaline 273-348 
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Table 2 
Systems containing C,H, with medium size components 




















2nd Component T [K] range P [bar] range ; LCVM 





Ay x 1000 AP% Ay x 1000 

















n -Heptane 235-500 ’ 43 
n-Octane 273-373 42 
n-Decane 278-478 5.2 
Methylcyclohexane 313-473 ; 5.4 
Propylcyclohexane 313-473 73 
Benzene 273-533 : 2.4 
Toluene 313-473 40 
Propylbenzene 313-473 , 9.7 
m-Xylene 313-473 ‘ 5.3 
Mesitylene 313-473 9.8 














Table 3 
Systems containing C,;Hg with medium size components 


























2nd Component T [K] range P [bar] range LCVM 





Ay x 1000 AP% Ay x 1000 























n-Butane 353-413 10 0.5 10 
n-Hexane 333-473 ’ 8 1.4 9 
n-Decane 278-S11 11 3.2 
Methylcyclohexane 313-473 6 3.8 
Propylcyclohexane 313-473 4.6 
Benzene 311-478 4.9 
Toluene 394-473 : 3.1 
Propylbenzene 394-473 4.8 
m-xylene 313-473 2.2 
Mesitylene 313-473 3.2 
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Table 4 
Prediction of the total liquid and vapor volumes (% devia- 
tions) for systems of CO, and n-alkanes 


LCVM 











2nd Component Ref. PHCT 





AV'% AVYY AVY 








Propane [20] 1.8 3.1 1.4 
n-Heptane [20] 3.0 5.4 3.1 
n-Octane [21] 4.0 10.0 10.3 
n-Decane [22] 1.2 9.1 an 








tion process involving petroleum cuts and natural 
products, is essential. Due to the complexity of 
the mixtures, which is very often in such cases, lack 
of experimental information is encountered. There- 
fore, in order to predict the phase behavior 
of such mixtures, fully predictive models are re- 
quired. The pure component parameters should be 
known, either experimentally, or by using correla- 
tions or other predictive methods. The binary 
interaction parameters should also be known. Such 
models are the LCVM model [1] and the PHCT 


Table 5 
System: CH,/nC,, [28] 








N T [K] range P [bar] range 














Table 6 
System: C,H,/nCj, [29] 





N T [K] range P [bar] range 
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equation of state [2]. A comparison of the ability 
of these two models to predict the phase behavior 
of such systems is of great interest. The fact that 
the development of these two models is derived 
from completely different concepts makes this 
comparison even more interesting. 

For the LCVM model, a predictive G* model 
(UNIFAC-type models) has been incorporated 
into mixing rules expressions for the attractive 
parameter of cubic equation of state (EoS) and it 
allows description of VLE in complex systems 
through a simple EoS. The pure component 
parameters are estimated from the critical proper- 
ties and the acentric factor (experimental or calcu- 
lated from correlations). The binary interaction 
parameters that are used, either already exist (e.g. 
UNIFAC interaction parameter table), or new 
parameters are estimated from binary systems. 

The PHCT EoS, considered in this work is based 
on the generalized van der Waals theory. The pure 
component properties are predicted from the 
molecular structure of the components of the 
mixture. The binary interaction parameters are 
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Table 7 
System: C,H,/nCoo [30] 





N T [K] range P [bar] range 
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Table 8 
System: C,H,/nC,, [31] 





N T [K] range P [bar] range 











Table 9 
System: C,H,/nCz,4 [32] 





N T [K] range P [bar] range 











estimated from mixtures and they represent 
families of components, since they are defined as 
segment—segment interactions. 

The common advantage of these two models is 
that the binary interaction parameters, which are 
essential for a correct prediction of the multicom- 


ponent phase behavior, can be calculated ‘a pri- 
ori’ by knowing the functional groups of the 
components present in the mixture. Consequently, 
these two models can be used as powerful predic- 
tive tools in the study of the behavior of such 
complex mixtures. 
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Table 10 
Systems of C,H, with large n-alkanes 





2nd Component T [K] range P [bar] range 





n-Hexadecane 
n-Eicosane 
n-Docosane 
n-Tetracosane 


320-360 8-114 
320-400 3-168 
320-360 2-93 

320-340 5-123 





* PHCT applies only in the far critical region (x,>0.2) 


2. Models 
2.1. The LCVM model 


Several equation of state/excess Gibbs energy 
(EoS/G®) models have been proposed for the 
prediction of vapor—liquid equilibria. LCVM, [1], 
MHV?2, [3] and PSRK, [4], [5] are three very 
successful models with wide range of applicability. 

The main characteristic of these EoS/G* mod- 
els is that the mixing rule for the attractive 
parameter « of the EoS is derived by setting the 
expression for G*, obtained from the EoS, equal 
to that from an existing G* model. This method 
can be applied at infinite pressure [6], or at zero 
pressure [7,8]. Boukouvalas et al., [1] have devel- 
oped a mixing rule that is a linear combination of 
the Vidal and Michelsen mixing rules. Of course, 
there is no need for specified reference pressure. 
The use of this mixing rule with the t-mPR EoS 
[9] and the Original UNIFAC [10] leads to the 
LCVM model studied here. 


2.1.1. The equation of state 

The EoS used is the translated and modified 
Peng—Robinson (t-mPR) EoS, [9]: 

RT 
V+ t—b 





Dw 


a 
~— (V+t) (V+t+b)+b V+t—b) 





The pure component parameter a is: 


2 2 


R 
a = 0.45724 I 


= en) 


4 

















For non-polar compounds (gases and alkanes), 
f(T,) 1s related to the acentric factor w: 


fT.) =[1l+4m(1 —-— /T,)}? (3) 
m = 0.384401 + 1.52276@ — 0.21380 2 
+ 0.0346 3? — 0.001976a 4 (4) 


Compared with the PR EoS using a single 
expression of m = f(@), similar to Eqs. (1) and (4) 
provides similar overall vapor pressure (P*) pre- 
dictions, but better ones at low P* levels. Addi- 
tionally, it gives very good predictions of 
saturated liquid densities, especially for com- 
pounds with large w values, through the introduc- 
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Fig. 1. Prediction of the phase behavior of CH,/nC,, system at 
320 K [28]; comparison between PHCT, LCVM and PR. 
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Fig. 2. Prediction of the phase behavior of C,H,/nC,, system 
at 320 K [30]; comparison between PHCT, LCVM and PR. 


tion of the translation. The expression for transla- 
tion t is given by Magoulas and Tassios [9]. The 
pure component parameter 5 is given by: 


RT. 
b = 0.0778 —— (5) 


Cc 


2.1.2. The mixing rules 
The mixing rule for the EoS parameter «a (= a/ 


bRT), is given by: 


RT 


A =) 


= Inat(Q+ Anu 
l1— A, b 
+ 4, Eain( > | (6) 


with 2=0.36, Ay = —0.623 and A,, = — 0.52. 
For parameter b, the conventional linear mixing 
rule is used: 


b = } x, b, (7) 


2.1.3. The G® model 
Original UNIFAC is used (in Eq. (6)), as it is a 


widely applicable G* model, which is being up- 

dated and extended regularly. The interaction 

parameters used in this work, were taken from: 

1. Hansen et al., [10], a recent revision of Origi- 
nal UNIFAC (temperature independent inter- 
action parameters, ‘¥;,=exp{ — A,/T}), and 

. Voutsas et al., [13] (linear temperature depen- 

dency, ,, = exp{ — [Ay + By: (7 — 298.15)]/ 
T)}). 

As demonstrated in previous publications 
[1,11—15] all models mentioned before present a 
comparable and very satisfactory behavior in the 
prediction of VLE in systems with components of 
similar size. The main application of EoS/G* 
models, however, is in the prediction of high 
pressure VLE, where the presence of gases with 
medium to high molecular weight compounds (i.e. 
significant size difference) represents the typical 
case. In this case, LCVM yields very good results, 
while the other models perform poorly. 

LCVM has been applied with considerable suc- 
cess to a variety of systems, such as: 

e CO,, CH,, C,H,, C3Hs, CsHy, N>, H2, HS 
with n-alkanes up to 44 carbon atoms and high 
pressures (up to 2000 bar). 

Other supercritical mixtures (e.g. gases with 

aromatics) 

Polar systems 

Multicomponent systems (ternary, gas conden- 

sates, oil systems, H.,S-rich sour natural gas, 

synthetic polydisperse mixtures) 

Solid—Fluid equilibria 


2.2. The PHCT EoS 


PHCT is an equation of state derived by Cot- 
terman et al., [16,17] based on the generalized van 
der Waals theory. The model has been developed 
in terms of the Helmholtz function considering an 
attractive and a repulsive contribution. The for- 
mer is made up of two terms, one for describing 
the high-density and one for the low-density re- 
gion, the latter being based on a generalized form 
of the Carnahan-—Starling equation. Dispersion 
and polar contributions to the free energy are 
considered in each part of the attractive term as 
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well as dipole and quadruple contributions. An 
interpolation function between the two limiting 
cases of a dense fluid and zero density fluid is 
included in the model based on empirical non-ad- 
justable parameters. This modification of PHC 
theory improves agreement with experimental sec- 
ond virial coefficients, vapor pressures, and satu- 
rated liquid densities. Since all molecular 
parameters have a well-defined physical signifi- 
cance, they can be reliably estimated for high- 
molecular-weight-fluids where experimental data 
are scarce. More important, separation into low- 
density and high-density contribution allows sepa- 
rate mixing rules for each density region; this 
flexibility in mixing rules significantly improves 
representation of mixture properties. Detailed ex- 
pressions about PHCT model are presented by 
Cotterman et al., [16,17]. 

The pure component parameters used, are pre- 
dicted, on the basis of the molecular structure, by 
means of the generalized correlations given by 
Gregorovicz et al., [18]. A fine-tuning of the 7* 
parameter may be performed if at least a single 
vapor pressure datum is available. 

The model includes at least three interaction 
parameters. One of them (that on the second 
virial term) becomes significant only at very low 
pressures, and, therefore, will be neglected in this 
investigation. The use of two unlike k,s is not 
justified by the different behavior in the two limit- 
ing regions for the systems investigated in this 
work, namely mixtures of hydrocarbons (non-po- 
lar) and supercritical gases. This reflects the sym- 
metric situation of the two infinite dilution 
regions, which are typical for non-polar systems. 

In conclusion, only one binary interaction 
parameter is used for a family of components, 
which accounts for the non-applicability of the 
random mixing in terms of segments making up 
the mixture (Fermeglia and Kikic, [2]). For polar 
systems at least two parameters are necessary. 

This model has been applied successfully to 
systems: 

e CO,, C,H, and C;H, with medium size com- 

pounds (up to 10 carbon atoms), [2] 

e CO,, CH,, C,H, with n-alkanes up to 44 car- 

bon atoms to the far-critical region, [17,19] 


3. Results 


The comparisons of these two models will be 
presented in two sections. In the first section, 
results for supercritical gases with medium size 
components will be presented. In the second one, 
results for supercritical gases with large n-alkanes 
will be presented. 


3.1. Mixtures with light gases and medium size 
components 


The supercritical gases are the CO,, C,H, and 
C,H,. The medium size components are alkanes 
(from 3 to 16 carbon atoms), methyl- and propyl- 
cyclohexane and aromatics (benzene, toluene, 
ethyl- and propylbenzene, m-xylene, mesitylene). 
The results that will be presented are VLE and 
volumetric predictions. Bubble point pressures 
predictions with LCVM and PHCT are presented 
in Tables 1—3. Volumetric predictions of these 
models are presented in Table 4. 

The results, presented in these tables, show that 
both LCVM, and PHCT with interaction parame- 
ters from Fermeglia and Kikic, [2], perform very 
well, even in their predictions of the volumetric 
behavior. In general terms, PHCT is slightly bet- 
ter in bubble point pressure predictions than 
LCVM, but two points must be stressed out: 

e LCVM interaction parameter have been ob- 
tained from correlation in the temperature 
range of 273-423 K. In some cases, it was 
necessary to extrapolate the use of these inter- 
action parameters to higher and lower temper- 
atures (from 234 up to 664 K). Although these 
results are not as good as the results obtained 
for the proper temperature range, they are also 
satisfactory and, therefore, the interaction 
parameters can be extrapolated for reasonable 
temperature range. 

All these systems have been used for the evalu- 

ation of the PHCT interaction parameters. On 

the other hand, for the evaluation of the 

LCVM interaction parameters only the gas/n- 

alkanes systems have been used. In the case of 

propane, the results with LCVM are straight 
predictions with the parameters suggested from 
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the Original UNIFAC revision by Hansen et al., 
[10]. 


3.2. Mixtures with light gases and large 
n-alkanes 


In this section, the two models are compared in 
terms of VLE calculations for systems containing 
gases and large n-alkane. The gases are the CHy,, 
C,H, and the n-alkanes are C,,¢, Coo, Co. and Cy,. 
These results are presented in the Tables 5—10 
and in the Fig. 1 and Fig. 2 

In this case, the interaction parameters used 
with the PHCT are the parameters proposed by 
Cotterman et al., [17] and the results are also 
presented in the paper by Alessi et al., [33] under 
the model name of PHCT3. The comparison of 
these results showed that LCVM performs very 
well to all asymmetries and all pressures (see also 
concluded Table 10). On the other hand, PHCT’s 
performance is not so straight. In the far critical 
region, the PHCT performs well, however, not 
better than LCVM. In the near critical region, the 
PHCT fails to give reliable results. This perfor- 
mance is well presented in Fig. 1 and Fig. 2. For 
comparison purposes, the performance of the PR 
EoS with conventional mixing rule (A;) is also 
presented in these figures. The interaction 
parameters used with PR EoS, were obtained by 
Abdoul et al., [34]. The behavior of PR EoS is 
very satisfactory and similar to the behavior of 
LCVM. 


4. Conclusions 


The performance of the LCVM model and the 
PHCT EoS is investigated in this work for binary 
systems of medium to large hydrocarbons with 
the light gases. The results indicated that the 
quality of their performance depends on the com- 
ponent size differences. 

For systems with components of medium size 
difference, both models can be applied with sig- 
nificant success. No problem seems to exist in the 
prediction of the VLE, as well as the prediction of 
the volumetric behavior of such systems. 

For systems with components of large size dif- 


ference, in the far critical region, the general 
picture is almost the same. In this region, LCVM 
almost always predicts the behavior of the systems 
slightly better than PHCT. However, in the near 
critical region, PHCT fails to provide reliable 
results, or any results at all. On the other hand, 
LCVM behaves satisfactorily. As it is mentioned 
by Alessi et al., [33], it is hard for PHCT EoS to 
describe the retrograde condensation of natural 
gas. On the other hand, LCVM performs such 
calculations successfully enough as it is presented 
by Spiliotis et al., [11], Boukouvalas et al., [15]. 


Nomenclature 


EoS attractive term (cohesion) 
parameter 

EoS couolume parameter 
Excess Gibbs free energy 
Characteristic temperature (en- 
ergetic pure component 
parameter) 


k,, Binary interaction parameter 


y 


Table symbols 
~ Vapor phase composition not 
available 
Liquid concentration of compo- 
nent 2[mol.] 
Average absolute percent error 
in liquid volumes 
AVY Average absolute percent error 
in vapor volumes 
AP% Average absolute error’ in bub- 
ble- or dew-point pressure 
Average absolute error in vapor 
phase composition 
N Number of data points per 
system 
Reference number 


Ay x 1000 
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Abstract 


Supercritical fluid extraction possesses several advantages over traditional liquid-solvent-based extraction methods, 
including improved selectivity, expeditiousness, automation and environmental safety. Fluid-phase equilibrium data 
are essential for studying the viability of supercritical extraction and the design of the extraction columns often used 
for this type of separation. Resveratrol has been correlated with serum lipid reduction and inhibition of platelet 
aggregation, and its cancer chemopreventive activity has recently been reported. The extraction of resveratrol using 
supercritical CO, is not possible; it only extracts traces of product, probably due to its polar nature. So, it is necessary 
to use a co-solvent. Ethanol is appropriate because it is a polar solvent permitted in the food industry. This work uses 
different mixtures of ethanol—CO,, ranging from 5 to 15% of ethanol, to optimise the amount of resveratrol 
extracted. The solubility of resveratrol in the supercritical solvent was measured at 313 K and pressures ranging from 
80 to 140 bar. In order to correlate the results obtained, two equations of state were successfully used (Peng— 
Robinson (PR) and Soave—Redlich- Kwong (SRK)). These equations show similar deviations, but PR shows a better 
correlation with experimental data. The solubility of resveratrol shows a maximum at an intermediate value of 
percentage of ethanol. It is roughly 7.5% of ethanol. So the optimal concentration of co-solvent was found. This 


solubility rises with pressure. © 2001 Elsevier Science B.V. All rights reserved. 
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parameter for cubic EOS (N m‘* mol ~ ”) 
parameter for cubic EOS (m? mol ~ ') 
binary interaction parameter in Eq. (3) 
number of data in Eq. (5) 
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universal gas constant (8.314 J mol~'!K~') 


absolute temperature (K) 
mole volume (m? mol ~') 
mole fraction 


Greek letters 
a acentric factor 
Ni binary interaction parameter in Eq. (4) 


Subscripts 

c critical property 

a component i, j respectively 
cal calculated results 

exp experimental data 








1. Introduction 


Supercritical fluids have been proposed for im- 
portant application in the field of extraction, anal- 
ysis, materials and reactions [1,2]. In fact, under 
proper conditions of pressure and temperature 
they have a liquid-like density and consequently 
quite good solvent properties. Supercritical fluid 
extraction presents several advantages over tradi- 
tional liquid-solvent-based extraction methods 
[3—5], including improved selectivity, expeditious- 
ness, automation and environmental safety. Selec- 
tivity can be significantly changed by varying both 
the pressure and temperature. The avoidance of 
organic solvents is a major goal in the isolation of 
natural products that may be commercialised as 
food additives. There are some companies in dif- 
ferent countries that are already working with 
carbon dioxide extraction processes, especially in 
the food and pharmaceutical industries where the 
quality of a product is no less important than 
considerations of cost. 

Phenolic compounds from plants are substances 
with the same metabolic origin; they are delivered 
from the shikimate pathway and phenylpropanoid 
metabolism [6]. The study of these compounds 
has gained special interest because of their phar- 
macological and, in some instances, cancer-pre- 


ventive properties. Resveratrol has _ been 
correlated with serum liquid reduction and inhibi- 
tion of platelet aggregation [7], and its cancer 
chemopreventive activity has recently been re- 
ported [8]. Phenolic compounds may be synthe- 
sised by’ plants for defence against 
microorganisms or strong UV radiation. For in- 
stance, ¢-resveratrol has been correlated with re- 
sistance to fungal infections [9,10]. 

t-Resvertatrol has been determined in wine, 
juices, berries and by-products of the grape-juice 
industry, including pomace, purees and seeds, by 
high performance liquid chromatography [11,12] 
and gas chromatography [13,14] in order to estab- 
lish the amount of resveratrol incorporated into 
the human diet and to evaluate these by-products 
as potential ingredients in food [12]. A combina- 
tion of geographical origin, variety, growing 
methods and winemaking procedures seems to 
affect the resveratrol concentration in wines and 
grape skin [10,11,14]. The content of t-resveratrol 
in berries steadily decreases during fruit ripening 
[13]. Grape skin is an important by-product of the 
wine industry. Now, it is used as animal food, but 
its possible use as a source of food ingredients in 
the human diet (as colourings, preservers or prod- 
ucts with pharmacological interest as resveratrol) 
is being investigated. 
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High-pressure vapour—liquid equilibrium data 
for the CO, + resveratrol system are required for 
rational equipment design. There are some tests of 
supercritical extraction of resveratrol in the litera- 
ture [15,16] that show the extraction of resveratrol 
using supercritical CO, was not possible. It only 
extracts traces of product. That is surely due to its 
polar nature. So, it is necessary to use a co-sol- 
vent. Ethanol is appropriate because it is a polar 
solvent permitted in the food industry. 

This work has as its objective the determination 
of solubilities of resveratrol in a mixture of etha- 
nol + CO,. This is a preliminary step in the study 
of the viability of supercritical extraction of 
resveratrol from grape skin in order to increase its 
presence in wines or for pharmacological use. 
Also, this information is necessary in order to 
optimise the extraction operating conditions. 


1.1. Models 


In order to correlate the results obtained, two 
equations of state were tried: the Peng—Robinson 
equation of state (PR EOS) [17], Eq. (1), and the 
Soave—Redlich—-Kwong (SRK EOS) [18], Eq. (2). 


a(T) 
v* + 2bv — b? () 








Trans-resveratrol 


Fig. 1. Structure of resveratrol. 


Table | 
Critical constants and acentric factor of pure components 





Substance T. (K) P.. (bar) 7) 





CO, 304.2° 73.7 0.2254 
Ethanol 513.92 61.4* 0.644* 
Resveratrol 893.8 51.4° 0.961° 





* Daubert and Danner [21]. 
> Estimated by the Ambrose method [22]. 
° Estimated by the Han and Peng method [23]. 


RT a(T) 





~v—b v(v+b) 


where v is molar volume; a (a measure of inter- 
molecular forces) and b (concerned with size of 
rigid spheres) are parameters of the equations, P 
is pressure, 7 is temperature and R is the univer- 
sal gas constant. 

The parameters were calculated by using the 
critical properties and the acentric factor of the 
pure components, and the following conventional 
mixing rules: 


(3) 
(4) 


where a,, a;, b; and b; are the parameters for pure 
components (function of critical properties and 
the acentric factor), and k,; and 4, are binary 
adjustable interaction parameters. 


2. Experimental section 


2.1. Chemicals 


Resveratrol (Fig. 1) (99 mass%, GC grade) was 
purchased from Aldrich Chemical Co., Inc. Etha- 
nol absolute (99.8 mass%, GC grade) supplied by 
Prolabo S.A. was used as modifier and solvent to 
collect the extract. The reagents were used with- 
out further purification after chromatography 
failed to show any significant impurities. High-pu- 
rity CO, (more than 99.9 vol% purity, SFC grade) 
was used as received. Calcined diatomaceous 
earth, extra pure (Scharlau Chimie S.A.), was 
utilised as a support for the products. The proper- 
ties of the pure components are listed in Table 1. 


2.2. Equipment and procedure 


Solubilities were measured by an SFX 3560 
Extractor with two Model 260D syringe pumps 
manufactured by ISCO (Lincoln, NE, USA). 
With the proper plumbing, the two-pump system 
can deliver modified supercritical fluid or a con- 
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Fig. 2. Schematic diagram of SFX 3560. Supercritical fluid extraction system. 


tinuous flow of supercritical fluid. The cylinder 
capacity of the pump is 266 ml and the maximum 
pressure 510 atm. The temperature can vary from 
313.2 to 423.2 K, and the supercritical fluid flow 
between 0.5 and 5 ml/min. The SFX 3560 sample 
reel holds up to 24 sample cartridges and collec- 
tion vials. The software supplied with the extrac- 
tor employs a dynamic method programming 
scheme which allows entry changes to methods 
while a previously programmed method is 
running. 

A schematic diagram of the apparatus is 
shown in Fig. 2. The sample cartridges with the 
product were placed in the cartridge reel, and the 
glass vials were taken up in the vial reel. Then 
the vial was raised, via a piston (vial lift), inside 
the collection system. Hence, the restrictor was 
inside the vial. After that, the cartridge lift raised 
the sample cartridge into the extraction chamber. 
Later, the vial was filled with the solvent, using 
the solvent pump, and simultaneously the extrac- 
tion chamber was pressured with CO, via the 
supply valve and preheat. After a static time 
(30s), the analyte valve was opened and the dy- 
namic time started (4—8 min). During this time, 
the supercritical fluid flowed through the product 
in the sample cartridge. CO, was dissolved in the 


solvent (vial) to trap the product. The wash valve 
allows washing of the tubes, using supercritical 
CO,, between the exit of the extraction chamber 
and vial. Finally, the vent valve was opened to 
decompress the extraction chamber. 

The extractor measured the amount of super- 
critical CO, and modifier ethanol used during the 
extraction. The flow of supercritical fluid used in 
all experiments was 2 ml/min. This is the suitable 
flow to achieve the equilibrium with saturation of 
supercritical phase. Ethanol was used to trap the 
resveratrol extracted. After that, the extract was 
transferred to 10 ml calibrated flasks and diluted 
with ethanol. In some cases, a higher dilution 
was necessary to achieve a concentration level 
adequate for UV spectrometry determination. 

For sample preparation, first 0.05 g of resvera- 
trol was dissolved in 10 ml of ethanol. Second, 
2ml of solution was added of diatomaceous 
earth (1 g) previously weighed in the sample car- 
tridge. Finally, the spiked solid was heated 
(60°C) for at least 2h in order to evaporate 
ethanol and precipitate resveratrol over the sur- 
face of the solid particles. This procedure is nec- 
essary in order to avoid channelling during 
extraction and the formation of masses of parti- 
cles. 
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An 8453 Hewlett-Packard diode array UV—VIS 
spectrophotometer was used for the determina- 
tions of resveratrol in the extract, after calibration 
with gravimetrically prepared standard solutions. 

Solubility is measured as the slope of resvera- 
trol moles extracted versus solvent moles (CO, + 
ethanol) at different extraction times. 


3. Results and discussion 


The extraction of resveratrol using supercritical 
CO, was tested. At extreme conditions of pressure 
and temperature it only extracts traces of product 
(<10~’ mole fraction), in accordance with Tena 
et al. [15]. This could be due to the polar nature 
of resveratrol. So, it is necessary to use a co-sol- 
vent. Ethanol is appropriate because it is a polar 
solvent and it is permitted in the food industry 
because it is harmless. Also, it is easily removed 
from the extract by evaporation at room 
temperature. 

The solubilities of resveratrol in the supercriti- 
cal phase were measured at 313 K, pressures rang- 
ing from 80 to 140 bar and for different mixtures 
of ethanol—CO,. These results are shown in Table 
pA 

This study used different mixtures of ethanol— 
CO,, ranging from 5 to 15% ethanol, in order to 
optimise the amount of resveratrol extracted with 
the minimal consumption of solvent. So, the ex- 
tract is cleaner and the cost of the solvent and 
separation is smaller. 


Table 2 


The composition of the supercritical phase was 
measured at 313 K. This is the minimum tempera- 
ture allowed for the equipment. This figure has 
been selected in order to avoid the opposed effect 
of temperature on solubility [1]. 

The lower limit of 80 bar was appropriate be- 
cause it was near of the CO, critical pressure, so 
the behaviour in conditions near the criticality of 
the system was studied. It was not possible to 
measure the solubility of resveratrol at pressures 
above 140 bar using this equipment because, at 
these conditions, the entire product contained in 
the sample cartridge was extracted. 

In order to correlate the results obtained, two 
equations of state were successfully used: the PR 
EOS [17] and the SRK EOS [18]. The parameters 
of the equations of state were determined using 
the solver function of Microsoft Excel and the PE 
program [19]. These parameters and the average 
relative deviation (ARD) are shown in Table 3. 
The parameters estimation was performed consid- 
ering all the experimental data by minimising the 
following objective function: 


OF = 5 ABSOen pe ‘ett) (5) 





where y is the molar fraction of solute (resvera- 
trol) in the supercritical phase. 

Generally, both the PR EOS and the SRK EOS 
show similar deviations at different conditions of 
pressure and temperature, but the PR EOS [17] 
shows a better correlation with experimental 
points. 


Solubilities (10°y3) of resveratrol expressed as mole fraction, in supercritical phase (CO, (1)+ethanol (2)) at 313 K and different 


percentages of modifier 

















P (bar) Resveratrol solubility at % modifier (ethanol (2)) of: 











5+0.01 7.5 +0.01 


10+ 0.01 








80+ 0.1 
90+ 0.1 
100 + 0.1 
110+0.1 
120+0.1 
130+0.1 
140+ 0.1 


4.64 + 0.07 
5.41 + 0.07 
6.33 + 0.07 
7.35 + 0.06 
8.47 + 0.02 
9.81+0.02 
11.24+0.01 


11.13 +0.02 
15.19 + 0.02 
18.63 + 0.02 
23.27 + 0.01 


9.31 + 0.03 
12.14+ 0.02 
15.08 + 0.02 
19.72 + 0.02 
27.01 + 0.02 


12.5+ 0.01 


15+ 0.01 




















5.69 + 0.07 
6.95 + 0.06 
10.20 + 0.03 
13.24 + 0.04 
19.07 + 0.02 


4.14+ 0.08 
5.10 + 0.07 
5.94 + 0.06 
7.59 + 0.06 
10.51 + 0.05 
13.66 + 0.03 
17.54 + 0.02 
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Table 3 


Parameters and deviations between experimental and calculated values of the system resveratrol (3)+ethanol (2)+CO, (1) for 


different equations of state at 313.15 K 





Model K1> 


N>3 ARD* 





PR EOS 0.075 
SRK EOS 0.019 


0.105 
0.116 





re ARD = (ry —Peall¥xn 


35 isonet raster eaenlmieahironrbaiadsnvia 








% ethanol 


Fig. 3. Solubility y of resveratrol (3) versus ethanol (2) fraction 
at: (MI) 80 bar, (A) 90 bar, (@) 100 bar; (@) 110 bar; (C2) 
120 bar, (A) 130 bar and (©) 140 bar. Modelling by PR EOS 
[17] shown by the line. 


As is shown in Fig. 3, the evolution of the 
solubility of resveratrol versus the modifier frac- 
tion shows a maximum at roughly 7.5% ethanol. 
This concentration allows the amount of solvent 
in the extract to be small; in this way the separa- 
tion process will be favoured. So the optimal 
concentration of co-solvent was founded. The ex- 
istence of a maximum for the evolution of the 
solubility with co-solvent concentration has al- 
ready been observed in the literature [20] in simi- 
lar cases. 

The solubility, as is shown in Table 2, rises with 
pressure. The extraction will be carried out at 
elevated pressure. The extraction pressure will be 
limited by operational costs and security cautions. 

In conclusion, the present study reveals the 
influence of process conditions on the solubility of 
resveratrol in supercritical CO,. It was necessary 
to use a co-solvent. Ethanol was an appropriate 
co-solvent and its optimal concentration was 
found to be about 7.5%. The solubility rises with 


pressure, so the extraction will be carried out at 
elevated pressure. Finally, the PR EOS [17] was 
applied successfully to the ternary system 
resveratrol—ethanol—CO,. 
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Abstract 


The use of biocatalysts in supercritical CO, (SCCO,) has received widespread attention in recent years. Biocatalysts 
have the advantage of substrate specificity under mild reaction conditions and SCCO, has several advantages over 
liquid solvents such as high solute diffusivities and low viscosity, which can accelerate mass transfer-limited enzymatic 
reactions. Concerning the enzymatic alcoholysis of vegetable oils, very little experimental data have been reported in 
the literature. Also, to the best of our knowledge, there is no comprehensive study available comparing conventional 
and supercritical media for this reaction. In this context, the main objective of this work is to compare enzymatic 
ethanolysis of palm kernel oil taking place in both SCCO, and in n-hexane as solvents. For this purpose, a Taguchi 
experimental design with two levels and four variables was adopted for each system to allow the investigation of the 
influence of process variables on the reaction conversion and on the enzyme activity. © 2001 Elsevier Science B.V. 


All rights reserved. 
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other ester (interesterification), resulting in a rear- 
rangement of the triglyceride fatty acid groups to 


Lipases have been extensively used in triglyce- produce a new triglyceride as a consequence of 
ride technology, mainly for the biotransformation the competitive hydrolysis and _ esterification 
of oils and fats. Among several important pro- reactions. 
cesses for lipid modification are the hydrolysis Esters obtained from alcohols and fatty acids 
reactions, synthesis of esters and transesterifica- have many important applications. For example, 
tion of these materials, catalyzed by lipases. In those from long chain acids (12—20 carbon atoms) 
these reactions the triglyceride reacts with a fatty and short chain alcohols (3—8 carbon atoms) have 


acid (acidolysis), an alcohol (alcoholysis) or an- been widely employed in food, cosmetics and 
pharmaceutical industries [1]. Natural esters such 


as those from jojoba oil, carnauba wax and whale 
oil have also been used. However, these oils are 


1. Introduction 
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expensive, and they are not usually available in 
large quantities. Therefore, it is desirable to de- 
velop methods for the production of such esters 
using cheaper and more plentiful raw materials 
[2]. Among several Brazilian raw materials of 
interest in the production of high-value-added 
products, palm fruit is one of the most prominent. 
It is commonly used for the production of edible 
oil and the kernel oil is used in the manufacture of 
soaps. Though direct application of alkyl esters is 
scarce, they can be used as intermediates in the 
oleochemical production. Note that this product 
can also be used as biodiesel [3]. 

Several researchers have reported an alternative 
method to produce esters through enzymatic reac- 
tions using lipases as catalysts [4-10]. Because 
biocatalysts have high specific activity and a low 
impact on the environment, they have become 
increasingly important for industry. For example, 
immobilized lipases are used as catalysts for reac- 
tions involving biomodification of triglycerides 
[11]. In this case, supercritical fluids (SCF), partic- 
ularly carbon dioxide, have currently received 
widespread attention as a possible medium for 
enzymatic reactions [12]. The main advantage of 
SCF over liquid solvents, such as n-hexane, is that 
the high diffusivity, low viscosity and low surface 
tension of SCF can speed up mass transfer-limited 
enzymatic reactions. 

Nevertheless, enzymatic ethanolysis of veg- 
etable oils has hardly been discussed in the litera- 
ture. Furthermore, no comprehensive study 
comparing conventional and supercritical media 
has been reported for this type of reaction. In this 
context, our main objective is to provide a com- 
parison between an organic liquid solvent and 
SCCO, for the enzymatic alcoholysis of palm 
kernel oil. For this purpose, a two-level Taguchi 
experimental design was adopted so as to investi- 
gate the influence of process variables, namely 
temperature (7), water and enzyme concentra- 
tions (denoted by [W] and [E], respectively) and 
oil/ethanol molar ratio (R) for n-hexane and tem- 
perature, water concentration, oil/ethanol molar 
ratio and pressure (P) for CO, system on the 
reaction conversion. The loss of enzyme activity 
for each experimental condition is also reported 
as well as the unusual effect of pressure on the 
initial reaction rate. 


2. Materials and methods 


2.1. Chemicals 


Palm kernel oil was used as purchased without 
any pre-treatment. Fatty acid composition of 
palm kernel oil was determined using a gas chro- 
matograph (HP 5890) with flame ionization detec- 
tor. The following instrumentation and conditions 
were used: H, as carrier gas, modified 
polyethylene glycol column (FFAP 2—25 m x 0.20 
mm i.d. x 0.30 um film), column temperature: 
180—210°C (2°C/min), injector temperature: 
250°C and detector temperature: 280°C. The fatty 
acid composition in palm kernel oil was deter- 
mined to be 9% caprylic acid, 47% lauric acid, 
15% myristic, 9% palmitic acid and 20% oleic 
acid, resulting in an average molecular weight of 
the oil as 701.9 g/gmol. Ethyl alcohol and n-hex- 
ane, both analytical grade, were used as substrate 
and organic solvent, respectively. Carbon dioxide 
with purity higher than 99.99%, was used as sol- 
vent in the high-pressure reactions. 


2.2. Enzymes 


Two commercial immobilized lipases were 
kindly supplied by NOVO Nordisk Bioindustrial, 
Brazil (Araucaria, PR): Rhizomucor  miehei 
(Lipozyme IM) immobilized on a macroporous 
anion exchange resin (0.15 U/g, 4% water and 
diameter in the range of 0.2—0.6 mm) and Can- 
dida antarctica (Novozym 435) immobilized on a 
macroporous anionic resin (0.12 U/g, 1.4% water 
and diameter in the range of 0.3—0.9 mm). The 
optimum activities for both enzymes are achieved 
at 40—60°C for Novozym 435 and 30—70°C for 
Lipozyme IM (Novo Nordisk, Technical Report 
A-05948, 1992). 


2.3. Analytical method 


The glycerol content evolved during enzymatic 
alcoholysis was determined using the method de- 
scribed by Cocks and van Rede [13]. The reaction 
conversion was then calculated by determining the 
glycerol concentration assuming a maximum glyc- 
erol yield at the end of the reaction of 12% of the 
mass of oil [14]. 
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Table | 


Process conversions obtained for the enzymatic alcoholysis of palm kernel oil in n-hexane 








Exp. Experimental conditions" 


Conversion (%) 











T CC) 


R Lipozyme IM Novozym 435 











40 
40 
70 
70 
40 
40 
70 
70 
a 12.5 


OMAAINHNDNBRWN — 


‘3 the 34.7 
1:10 52.1 42.1 
1:3 28.9 54.5 
1:10 S757 53.4 
1:3 62.2 41.5 
1:10 34.2 58.3 
1:3 29.8 35.2 
1:10 31.9 29.2 
1:6.5 32.9 ES Be 








*[E], enzyme concentration, wt.% of substrates; [W], water concentration, wt.% of substrates; R, oil/ethanol molar ratio. 


2.4. Apparatus and experimental procedure 


2.4.1. Conventional solvent 

The experiments were performed in stoppered 
125-ml Erlenmeyer flasks. Lipase at pre-estab- 
lished concentrations was added to the mixture of 
oil (1 g)—ethanol and n-hexane (40 ml) and the 
flasks were agitated at 200 rev./min for 6 h in a 
controlled-temperature shaker. The organic sol- 
vent was chosen based on preliminary experi- 
ments with n-hexane and petroleum ether, which 
presented practically the same reaction conver- 
sion. The former solvent was then used as reac- 
tion medium due to its higher normal boiling 
point thus facilitating the experimental execution. 
The agitation level and the reaction time are 
suitable values obtained after a number of experi- 
mental tests. Indeed, an enhancement of these two 
operational variables meant negligible changes in 
observable conversions. A Taguchi experimental 
planning with two levels and four variables (tem- 
perature, water and enzyme concentrations and 
the oil/ethanol molar ratio) was adopted. The 
experimental planning, covering the variable 
ranges commonly employed in the literature [15], 
is presented in Table 1, columns 1 to 5. The 
experiments were accomplished randomly, and 
duplicate runs were carried out for each experi- 
mental condition resulting in a reproducibility 
better than 5%. The process conversion was then 
modeled by an empirical expression. 


2.4.2. SCCO, 

The experimental set-up, as shown in Fig. 1, 
consists basically of a CO, cylinder, a high pres- 
sure pump (ISCO 260D), a liquid pump (Thermo 
Separation Products, model constaMetric 4100) 
and a 300-ml reactor (Autoclave Engineers, model 
300BG) equipped with a mechanical stirrer (Mag- 
neDrive II), a heating mantle, an internal cooling 
loop, capable of keeping the temperature constant 
within 0.1°C, a pressure transducer and a sam- 
pling tube. Palm kernel oil at pre-established con- 
centrations, enzyme (5 wt.% of substrates) and 
water (0-10 wt.% of substrates), were placed in 
the reactor. Afterwards, the reactor was closed, 
flushed and pressurized with CO,. Before reaching 


















































Fig. 1. Schematic diagram of the experimental apparatus: A, 
CO, cylinder; B, syringe pump; C, reactor; D, temperature 
controller; E, pressure transducer; F, liquid pump; G, cold 
trap; H, check valve; I, micrometering valve. 
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Table 2 


Process conversions obtained for the enzymatic alcoholysis of palm kernel oil in SCCO, 








Exp. Experimental conditions 


Conversion (%) 








Cu P (bar) 


Lipozyme IM 


Novozym 435 








40 73 
40 73 
70 
70 
40 
40 
70 
70 
55 


OSCOAAIHD MN HWW — 


7.7 54.5 
16.7 63.2 
11.9 24.7 
14.2 35.4 
11.7 14.1 
25.4 21.3 
15.5 12.1 
18.4 42.9 
26.4 20.4 








the desired pressure, a pre-established amount of 
ethanol was fed into the reactor through the use 
of the metering liquid pump. The substrates and 
enzyme were then continuously mixed with an 
agitation level of 1000 rev./min in order to 
provide a proper homogenization of the reacting 
mixture. Moreover, it was experimentally ob- 
served that the degree of agitation did not affect 
sensibly the enzyme activity. The amount of al- 
lowable reactants were determined based on the 
phase equilibrium data from Bharath et al. [17] by 
performing simple approximate calculations so as 
to increase the liquid phase volume thus avoiding 
unnecessary ethanol charges and to permit reli- 
able samplings — these data are reported else- 
where [16]. 

To assess the influence of process variables on 
the reaction conversion, the reaction time was set 
at 4 h. Typically, samples of approximately 2 ml 
were withdrawn after the reaction completion, or 
at the beginning of the reaction, and then ana- 
lyzed using the previously described method. A 
Taguchi experimental design with two levels and 
four variables (temperature, water and enzyme 
concentrations and oil/ethanol molar ratio) was 
adopted. The experimental design along with the 
variable ranges is presented in Table 2, columns 
1—5. Duplicate runs were made for all experimen- 
tal conditions and the variations around the aver- 
age values were typically about 3%. The process 
conversion was then modeled by empirical 
expressions. 








3. Results and discussion 


The experimental results using n-hexane and 
SCCO, are presented in Tables | and 2, respec- 
tively. One can see from these tables that the 
highest conversion (63.2%) in SCCO, was ob- 
tained using Novozym 435 as catalyst while in 
n-hexane Lipozyme IM provided the best achieve- 
ment (77.5%). These results might be very impor- 
tant from an economic standpoint since a 
by-product can be used to produce high-value- 
added substances at near room temperature. 
Moreover, besides the good conversion obtained 
in SCCO,, it might be relevant to call attention to 
the much lower expense of solvent in supercritical 
medium, 2 ml CO,/g oil, compared with conven- 
tional solvent, 40 ml n-hexane/g oil. 

From the results obtained in the experimental 
design, a standard statistical modeling technique, 
using the Statistica® software, was employed in 
constructing empirical models in order to evaluate 
the effect of process variables on the reaction 
conversion. The empirical models were built as- 
suming that all variable interactions were signifi- 
cant, estimating the parameters related to each 
variable interaction and main variable effects and 
discarding the meaningless parameters by using 
the Student’s t-test considering a confidence level 
of 95%. The objective of using the Student’s t-test 
is to evaluate if the parameters were significantly 
different from zero. This test takes into account 
the S.D. of each parameter according to a well 
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known procedure available in many textbooks 
[18]. The parameters were estimated through the 
maximum likelihood method [19]. With regard to 
the parameter analysis, it is important to mention 
that a parameter with a negative value implies a 
negative effect of the variable on the process. 


3.1. Conventional medium 


It can be noticed from Table 3 that for the 
system containing Novozym 435, temperature, 
water and enzyme concentrations and also the 
cross interaction between temperature and oil/eth- 
anol molar ratio significantly affect the process 
conversion. Though the optimum temperature for 
this enzyme is about 70°C, a negative effect is 
observed on the conversion. The added water also 
presented a negative effect, corroborating that 
water excess may change the reaction equilibrium, 
diminishing the formation of esters. As expected, 
the enzyme concentration, within the range stud- 
ied, presented a positive effect on the conversion. 
Besides the isolated effects of these variables, 
cross interactions should also be taken into ac- 
count, since cross interaction between tempera- 
ture and oil/ethanol molar ratio had a negative 
effect on the conversion. It can be concluded that 
an excess of alcohol at higher temperatures leads 
to a decrease in conversion. 

With regard to the system containing Lipozyme 
IM, temperature and oil/ethanol molar ratio have 
the most pronounced effects. Temperature pre- 


Table 3 


sented a negative effect on the reaction conver- 
sion, the lowest value (40°C) being the condition 
in which the enzyme presents the best perfor- 
mance, according to the manufacturer. In this 
case, the enzyme concentration showed a negative 
effect on the reaction conversion. Table 3 also 
shows for this system that an excess of ethanol 
may inhibit the enzymatic reaction. It is also 
worth noticing the influence of the cross interac- 
tions on conversion, namely temperature—oil/eth- 
anol molar ratio and temperature—enzyme 
concentration. The cross interactions are very 1m- 
portant since from a practical standpoint the 
analysis of only isolated (main) effects may lead 
to unexpected/undesirable results. 


3.2. SCCO; 


Based on the results obtained in conventional 
solvent, the first system investigated in SCCO, 
was the alcoholysis of palm kernel oil catalyzed 
by Lipozyme IM. However, by comparison of 
Tables 1 and 2, the conversions obtained in this 
case were, for all experimental conditions, lower 
than those obtained in a conventional solvent. 
For this system, it can be observed from Table 4 
that the oil/ethanol molar ratio was the variable 
that more positively affected the process conver- 
sion. For the conventional medium, this was an 
expected result, since in supercritical medium, the 
ethanol is more susceptible to migrate to the gas 
phase. The cross interactions between tempera- 


Regression results (7 = d)+a,T+4,[E]+a,R+a,[W]+a;7T[E]+a,7TR) for the system containing Novozym 435 and Lipozyme IM in 


n-hexane* 





Novozym 435 


Lipozyme IM 











Effect Parameter 


Effect Parameter 





(do) 41.03 

(a,) —3.15 
[E] (a>) 6.01 

(a3) —0.09 
[W] (a4) —2.66 
T[E] (as) 0.38 
TR (dg) —6.31 


(do) 43.89 
(a,) — 12.05 
(a>) —3.23 
(a3) —5.07 
(a4) 1.39 
(as) 4.69 
(a) 8.09 











*, calculated conversion; ¢, parameter S.D. 
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Table 4 


Regression results for the system containing Novozym 435 and Lipozyme IM in SCCO, 





Lipozyme IM 
N =d)+a,P+a,R+a,[W]+a,TP+a;TR+a,PR+a,T° 


Novozym 435 
N =ayt+a,T+a,P+a,R+a,[W]+a;TP+a,TR+a,T* 





Effect Parameter 


Effect Parameter 





(dy) 24.27 
(a,) 0.61 
(a) 3.48 
(a;) —0.66 
(a4) —2.57 
(as) —2.19 
(as) 0.51 
(az) —9.1 


(do) 20.12 
(a,) —2.39 
(a5) —11.29 
(a3) 4.83 
(a4) —0.72 
(as) 8.76 
(a,) 1.80 
(az) 14.83 





Table 5 


Observed activity losses for the enzymatic alcoholysis of palm kernel oil in SCCO, 





Exp. Experimental conditions 


Loss of enzyme activity (%) 





T (°C) P (bar) 


Lipozyme IM Novozym 435 





40 73 
40 73 
70 200 
70 200 
40 200 
40 200 
70 73 
70 73 
55 


11.88 23.05 
25.05 4.9] 
35.81 41.23 
58.39 49.00 
38.35 38.74 
31.66 44.31 
41.74 35.31 
20.06 28.86 
74.26 45.61 





ture—pressure and temperature—oil/ethanol molar 
ratio also had a significant effect on the reaction 
conversion, both of them showing a negative ef- 
fect. From this comparative study, it becomes 
clear that enzymes may exhibit different behaviors 
depending on the solvent medium and thereby 
affecting the process conversion. 

In a second step, the same procedure was 
adopted using Novozym 435 as catalyst. After 
carrying out all the experiments, an empirical 
model was built for this system. It can be noted 
from Table 4 that the pressure was the variable 
that more strongly affected the process conver- 
sion, revealing that this variable has a negative 
effect on the conversion for this enzyme, probably 
causing enzyme activity loss. The oil/ethanol mo- 
lar ratio also had a significant effect, similar to the 


Lipozyme IM. Temperature also presented a sig- 
nificant negative effect and hence an increase in 
this variable may lead to inactivation of the en- 
zyme. The interaction between temperature and 
pressure also had an important effect on the 
conversion. Regarding the reaction conversion, as 
noted for n-hexane, different values have been 
recorded when using Lipozyme IM and Novozym 
435 in SCCO.,. 

While for the conventional medium no loss of 
enzyme activity was experimentally observed, the 
same was no longer valid for the SCCO, system. 
Steinberger and Marr [20] have recently pointed 
out that the stability of an enzyme in SCCO, 
depends on both its tertiary structure and on 
several parameters during exposure to high-pres- 
sure CO,. They contend that high temperature, 
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the water content in CO, and pressurization/de- 
pressurization steps can cause enzyme inactiva- 
tion. In fact, as presented in Table 5, in most of 
the experimental conditions considerable losses of 
enzyme activity have been observed. An attempt 
at explaining this is very difficult taking into 
account a complex reaction medium where many 
parameters are simultaneously involved. Thus, it 
seems difficult in this case to ascertain a relation- 
ship between enzyme activity and _ process 
conversions. 

Another interesting aspect is the effect of pres- 
sure on the enzyme-catalyzed reactions in super- 
critical fluids. It has been observed by Cernia et 
al. [21] that the initial reaction rate generally 
exhibits an unusual behavior in the range of 6—25 
MPa. This curious phenomenon was explored in 
this work using Novozym 435 as catalyst at the 
optimized values of process variables (T = 40°C, 
[W] = 0 and R = 1:9.8). As shown in Fig. 2, in the 
subcritical region a rise in pressure leads to a 
sharp increase in the reaction rate and then a 
gradual decrease at higher pressures. The rate of 
glycerol production rapidly increases from 13.09 
(60 bar) and 16.08 (73 bar) up to 28.15 at 90 bar 
and then dramatically decreases to 6.66 at 200 
bar. The fact that the initial rate of enzyme-cata- 
lyzed reactions is SCF is maximized near the 
solvent critical region has just been corroborated 
in the literature [21]. Though the role of SCCO, 


147 


molecules during enzymatic reactions is far from 
understood, it has been argued that the pressure 
could induce modifications of both structure and 
function of enzyme leading to a change in the 


enzyme activity [21]. 


4. Conclusions 


In this work, Lipozyme IM and Novozym 435 
were used as catalysts in the ethanolysis of palm 
kernel oil in n-hexane and SCCO.,. The use of an 
experimental design for the production of esters 
from enzymatic reactions of vegetable oils proved 
to be a rational basis to assess the influence of 
process variables on the conversion. From this 
study one can verify that depending on the en- 
zyme and solvent used, different results can be 
obtained. The main interest of this work lies in 
the fact that it is the first study comparing con- 
ventional and supercritical media for the alcohol- 
ysis of vegetable oils. From these results an 
investigation concerning the kinetics of the reac- 
tion in conventional and supercritical media can 
be carried out. Despite enzyme activity losses, 
SCCO, seems to be an appropriate solvent for 
enzymatic alcoholysis, very similar to an organic 
liquid solvent with regard to process conversion. 
Moreover, an advantage of supercritical fluids, 
considering the complete process, is that after the 
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Initial rate (“glycerol/min.) 














120 


150 


Pressure (bar) 


Fig. 2. Effect of pressure 
R= 1:9.8. 


on the initial reaction rate of palm kernel oil catalyzed by Novozym 435 at T= 40°C, [W]=0 and 
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system depressurization, a mixture of products 
and non-reacted substrates may be obtained with- 


out traces of solvent. 
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Abstract 


Extraction of egg yolk oil (EYO) from egg yolk powder (EYP) with supercritical CO, was performed on a 
laboratory apparatus. Solubility of EYO in supercritical CO, was measured. The external diffusion and the 
equilibrium between solid and fluid phases were experimentally found to be the controlling steps during the extraction 
process. Based on this mechanism, a mathematical model for this extraction process was developed. The model 
parameters, adsorption equilibrium constant (k,), EYO concentration in solid controlling the transition in the 
equilibrium (x,) and the overall volumetric mass transfer coefficient (K,a,), were obtained by simulation. The 
simulation results indicated that x, is 0.56, K,a, is directly proportional to CO, flow rate with an exponent of 0.548, 
and the adsorption heat of EYO is 6-9 kJ mol~'. The model was verified by concentration profiles of solid. The 
extraction process of EYO with supercritical CO, was conducted on a pilot plant and the developed model could 
predict satisfactorily the process. © 2001 Elsevier Science B.V. All rights reserved. 


Keywords: Extraction; Supercritical CO,; Egg yolk oil; Mathematical model 





[1]. Therefore it is important to remove choles- 
terol before obtaining high purity phospholipid. 
Traditionally phospholipid is extracted from EYP 
with organic solvent, such as acetone, ethanol and 


1. Introduction 


Egg yolk powder (EYP) contains abundant 
phospholipid, which is of high value in food, 


medicine and cosmetics. It also contains choles- 
terol, glyceride, phospholipid, protein and mois- 
ture. The high level of cholesterol in EYP is the 
major problem that limits the EYP application 
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E-mail address: wzwu@ms.sxicc.ac.cn (W. Wu). 


hexane [2]. However, this method has its in- 
evitable disadvantages such as organic solvent 
residues in products and degeneration of active 
components by the required high temperature. 
Recently, supercritical carbon dioxide (SC-CO,) 
extraction is applied in separation of EYP [3-6]. 
In this process, cholesterol, glyceride and mois- 
ture, called egg yolk oil (EYO), are dissolved in 
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SC-CO,, and phospholipid and protein, called 
‘matrix’, leave in the solid phase. Using this tech- 
nique, Froning et al. [3] and Hung et al. [4] 
reported cholesterol removal efficiency of about 
70% and Wu et al. [5] achieved a complete separa- 
tion under optimal separation conditions of 55°C 
and 28-36 MPa. The merits of this method are 
not only avoiding disadvantages mentioned in 
organic solvent extraction, but also not impairing 
the functional properties and taste of EYP [3,4]. 
Hence this technology has potential application in 
commercial plant. Up to now, no model has been 
built for SC-CO, extraction of EYO from EYP. 
Therefore it is necessary to study the mass trans- 
fer mechanism and to build a mathematical model 
for this process. 

A few models had been established for super- 
critical fluid extraction (SFE) of other materials. 
For building models, different mass transfer 
mechanisms had been hypothesized and most 
models proposed were built on different mass 
balance integration. Lee et al. [7], Bulley et al. [8] 
and Fattori et al. [9] assumed that mass transfer 
resistance occurred only in fluid phase when they 
set up models of SC-CO, extraction of canola 
seed oil. Reverchon et al. [10,11] found that the 
controlling step was the internal diffusion in solid 
phase during extraction of herbaceous essence and 
sage oil. Roy et al. [12] made a distinction in solid 
phase between broken cells and intact cells, which 
were considered with different mass transfer resis- 
tance: when extracting oil in broken cells the 
resistance occurred in fluid phase and when in 
intact cells it did in solid phase. An extraction 
model with a shock transition concentration and 
variable mass transfer resistance was proposed for 
fractionation of flower concrete using SC-CO, 
[13]. Roy et al. [14] adopted the shrinking core 
model for extraction of ginger oil. Desorption 
equilibrium was considered as the limiting step 
during SC-CO, extraction of essential oils [15]. 
Based on above reviews, different models were 
built for SFE from different materials. 

The objectives of this work are: (1) to study the 
mass transfer mechanism of SC-CO, extraction of 
EYO from EYP; (2) to build a mathematical 
model for this process; and (3) to verify the model 
and to apply it to a pilot plant. 


2. Experimental 


2.1. Materials 


EYP, obtained from Beijing International Egg 
Products, contains EYO of 43.3 wt.%. EYP was 
stored in a refrigerator and was extruded into 
three grades of bars with diameters of 1.2, 2.0, 
and 3.0 mm, respectively, before fed into extrac- 
tor. CO,, bought from Shanxi fertilizer factory, 
China, has a purity of more than 99%. 


2.2. Experimental apparatus 


The experimental apparatus is shown in Fig. 1, 
which mainly consists of an extractor (E) and two 
separators (SI, S2). The extractor is 4.4 x 10~* m 
in diameter and 1.2 x 10~? m? in volume. Each of 
the two separators in series, with a volume of 
0.5 x 10~3m?, has a conical bottom connected 
with a regulating valve. A high-pressure plunger 
pump (PU) (Benxi water pump factory, Benxi 
China) is equipped for operating at pressure up to 
40 MPa and at CO, delivering flow up to 15 kg 
h~'. The flow rate of CO, was controlled by 
adjusting the pump speed and monitored by an 
electronic high-pressure rotameter (R) (Beijing 
forth automation meter factory, Beijing, China). 
Pressures in extractor and separators were con- 
trolled by means of manually adjusted microme- 
tering valves (Funing Gauge and Valve Factory, 
Funing, China). Temperatures in heat exchangers, 





Fig. 1. Schematic diagram of the extraction of EYO from 
(EYP) with SC-CO,. C, Cooler; CY, CO, cylinder; E, extrac- 
tor; F, filter; H, heater; PU, high pressure pump; ®) , pressure 
gauge; R, rotameter; Sl, separator 1; S2, separator 2; ® , 
thermocouple; TK, CO, tank. — — — heated water circle; 
cooled water circle; —-—- CO, circle. 
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extractor, separators and tank were controlled by 
water bathes that maintained the operating tem- 
perature within +0.5°C. Extractor temperature 
was monitored continuously with a thermocouple 
in extractor exit. The pressures, monitored at 
pump exit and at extractor exit, did not give a 
meaningful difference. Consequently the pressure 
monitored at extractor exit was considered as the 
operating pressure. 


2.3. Experimental procedure 


A desired sample was fed in extractor after 
some glass beads were put at the bottom of the 
extractor for uniform distribution of CO, flow. 
Air in the extractor was removed by CO,. Then 
the desired temperatures of extractor, separators, 
CO, tank and heat exchangers were maintained 
by water bathes, respectively. Liquid CO, from 
the CO, tank was pumped into the extractor 
through the heat exchanger. When reaching the 
desired value, the extractor pressure was con- 
trolled by adjusting the micrometering valve at 
the separator entrance and CO, was passed 
through the extractor column to start the extrac- 
tion of EYO from EYP. CO, bearing EYO ex- 
panded through the valve and lost its solution 
power to EYO. Then EYO extracted by CO, was 
precipitated at the separators and discharge of the 
separators was performed at fixed time intervals. 
CO, free EYO flowed into the heat exchanger 
through a rotameter, which measured the CO, 
flow rate, and became liquid CO, and returned 
into the liquid CO, tank. 


2.4. Analysis 


Content of EYO in EYP was analyzed using 
the method in reference [5]. The determination of 
cholesterol content in EYO follows a method as 
described in the literature [16]. The content of 
moisture was measured by a vacuum oven drying 
procedure (16.002, AOAC, 1984). The composi- 
tion of fatty acids in EYO was analyzed using a 
gas chromatograph (GC-14B Shimadzu, Japan) 
equipped with a flame ionization detector. The 
column is | m long, 3 mm internal diameter, 
carried polyethylene glycol succinate supported 


on white diatomaceous earth. Saponification and 
methylation of samples were based on method 
969.33, AOAC (1990). The methylated samples 
were analyzed under the column temperature of 
180°C, the injection temperature of 250°C and the 
detector temperature of 250°C. The fatty acid 
peaks were identified by comparison with the 
retention times of standard fatty acids assayed 
under identical conditions. 


3. Mathematical model 


3.1. Extraction mechanism 


Since the employed CO, flow rate is very low 
(<5~x10~7 ms~'), the EYP particles in extrac- 
tor can be considered to be stationary, and the 
extraction process is similar to a fixed-bed pro- 
cess. When CO, flows in the void of EYP particles 
in the extractor, EYO molecule is dissolved in 
SC-CO,. The extraction of EYO from EYP is an 
unsteady process. The extraction mechanism for a 
single particle is described by the following steps 
in series: (1) EYO molecules are desorbed from 
matrix, and dissolved in SC-CO,, which is called 
dissolution equilibrium process; (2) EYO 
molecules dissolved in SC-CO, diffuse through 
the solid holes to external surface, which is called 
internal diffusion; (3) EYO molecules diffuse from 
external surface through the SC-CO, film to fluid 
body, which is called external diffusion. 


3.2. General hypotheses 


To simplify mathematical model, following hy- 
potheses are needed: (a) temperatures and pres- 
sures at all parts of the extractor are uniform; (b) 
axial dispersion of EYO is negligible; (c) EYO is 
uniformly distributed in EYP particle; (d) the 
volume fraction of extraction bed, @, is not influ- 
enced by the weight loss of EYP during the 
extraction; (e) although several components are 
involved in the extract, their mass transfer behav- 
iors are supposed to be similar and to be de- 
scribed as a single pseudocomponent, which is 
partly based on the experimental results that the 
different components in the extract are not sepa- 
rated during the extraction (see Table 1). 





152 W. Wu, Y. Hou / J. of Supercritical Fluids 19 (2001) 149-159 


Table | 
Composition of egg yolk oil (EYO) extracted during different extraction periods* 





Sample” Composition of fatty acid in EYO/wt.% Cholesterol/wt.% Moisture/wt.% 








Ci). CURE (AiG. Cig: Cie  <iet: Cite: Crs 





0.50 0.20 25.26 5.02 6.81 45.45 15.84 0.92 6.99 
0.27 0.09 24.08 3.43 9.40 48.02 13.87 0.84 5.03 


0.41 0.14 24.63 4.54 


7.26 47.57 14.43 0.92 6.59 








“ Extraction conditions: 28 MPa, 55°C, 7.0 kgCO, h~', 0.250 kg of EYP with a particle size of 1.2 mm. 
> §-1, sample collected at first 30 min; S-2, at last 30 min; S-3, average sample. 


3.3. Mass balance equations 


According to the above hypotheses, the differ- 
ential balance of EYO in the fluid phase along the 
extraction bed is 

Oy oy 


i aoe a Se L 
pe a, + P U ah J(x, y) (1) 


and that in the solid phase is 
(2) 


where ? is the fluid density, which is supposed not 
to be influenced by the presence of EYO; @, is 
the density of the matrix (free EYO); y is the 
concentration of EYO in SC-CO,; ¢ is the con- 
sumed time; U is superficial velocity of CO, 
through the extractor; / is bed length; x is the 
concentration of EYO in EYP; J is overall volu- 
metric mass transfer rate. 

The initial condition on x and y and the 
boundary conditions on y for finding the solution 
of Eqs. (1) and (2) are following: 


i x(h, t=0)=x, 
< y(h, t=0)=0 
| y(t, h =0) =0 





3.4. Overall volumetric mass transfer rate, J 


The expression of J is experimentally deter- 
mined, but it is strongly dependent on the struc- 
ture of matrix and the behavior of EYO transfer 
in the solid and fluid. Therefor, it is first to 
understand the dependence. 


The effect of the external mass transfer on the 
extraction was studied via some experiments of 
fixing contact time (fixing the ratio of EYP 
amount to flow rate of CO,). The results is shown 
in Fig. 2, which shows that the rate of EYO 
extraction increases with increasing flow rate of 
CO,, and then suggests that the external mass 
transfer of EYO should be considered when J 
expression is built. 

Fig. 3 shows the effect of particle size on the 
extraction rate of EYO. There is no obvious 
change in the extraction rate of EYO at different 
particle sizes, which suggests that the internal 








—a— CF 11.98kgh 
—1— QE 301kgh 




















10 150 20 20 
Time/min 


Fig. 2. Effect of flow rate of CO, on extraction under the 
fixing contact time. Conditions: temperature, 45.0@; pressure, 
32.0 Mpa; particle size, 1.2 mm in diameter; CO, flow rate, 
3.01 kgh~! and 11.98 kgh~'! for 0.125 kg and 0.500 kg EYP, 


respectively. 
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Fig. 3. Effects of particle size of EYP on extraction. Condi- 
tions: temperature, 45.0@; pressure, 32.0 Mpa; CO, flow rate, 
about 7.0 kgh~'; EYP, 0.250 kg. 


mass transfer of EYO have little influence on the 
extraction process. 

In order to characterize the equilibrium of EYO 
dissolution between the two phases, the experi- 
ments at different CO, flow rates were performed 
and the results are shown in Fig. 4. The four sets 
of data exhibit no clear difference, which suggests 
that the amount of extracted EYO is constant 
when the amount of fed CO, is given and extrac- 
tion rate is dependent on the CO, flow rate. It 
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Fig. 4. Effect of CO, flow rate on extraction. Conditions: 
temperature, 45.0@; pressure, 32.0 Mpa; particle size, 1.2 mm 
in diameter; EYP, 0.250 kg. 


means that there is a relationship of concentra- 
tions between the two phases [11]. 

To determine the relationship, it is supposed 
that there are two forms of EYO molecular in 
EYP: one is strongly adsorbed on the matrix 
molecular, the other is weakly adsorbed on EYO 
molecular, and the extraction of the latter is prior 
to the former. The concentration of strongly ad- 
sorbed EYO is expressed as x,, then the concen- 
tration of the weakly adsorbed EYO is (x )—x,), 
where X, is the initial concentration of EYO in 
EYP. When the concentration of EYO in EYP, x, 
is greater than x,, the fluid concentrations at the 
contact face of solid, y*, reach the solubility of 


EYO, y,. That is 
yt=y,. (4) 


when x is less than or equal to x,, to evaluated the 
equilibrium relationship between y* and x, a lin- 
ear dependence on solid concentration, x, is 


adopted, 
(5) 


where k,, is the adsorption constant of EYO be- 
tween the solid and the fluid phases at equilibrium. 
Based on the description above, J can be ex- 


pressed as follows 


Kia,p(y,—y) at 72> x, 
a _ a Pp r 
thane =e Bagi —y)at xsx, 


(6) 


where a, is the specific surface area of EYP 
particle, K; is the external mass transfer coefficient. 
Substitution of (6) into (1) and (2) leads to a 
partial differential equation system: 


os Vr 
x/k, at X SX, 


After treated dimensionlessly, the equation sys- 
tem can be numerically solved with finite element 
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Table 2 
Solubilities of egg yolk oil (EYO) in SC-CO, (g of EYO/kg of 


CO.) 











Pressure (MPa) Temperature (°C) 








35.0 45.0 55.0 65.0 75.0 











pe 5.7 5.6 4.7 3.8 
6.9 6.9 7.4 6.8 6.6 
8.6 9.6 10.0 9.5 9.4 
9.9 11.4 13.1 13.2 














method using the boundary and initial conditions 
(3). 

As a result, two function, y= y(t,4) and x= 
x(t, h), can be obtained by solving Eq. (7) numer- 
ically. Therefore it is possible to calculate the 
EYO recovery, E(t): 


E(t)= ZA | y(t, A) dt, (8) 


where Q is the CO, flow rate and W is the amount 
of EYP fed in extractor (kg). X, is the concentra- 
tion of EYO in raw material. Thus E(t) is the 
ratio of the amount of EYO extracted to the 
amount of EYO in raw material. 


3.5. Determination of parameters for 
mathematical model 


Experimental extraction curves for measuring 
solubility of EYO in SC-CO, can be divided into 
two parts, named fast extraction and slow extrac- 
tion periods. The former process is controlled by 
the solubility of EYO in SC-CO,. The latter pro- 
cess maybe controlled by mass transfer or desorp- 
tion velocity. Solubility, y,, was evaluated from 
the slope of the extraction curve in the initial 
period [17—19]. The solubility of EYO in SC-CO, 
was measured with the same apparatus. The re- 
sults show in Table 2, which was carried out at a 
flow rate of about 4.0 kg h~!' and at various 
pressures and temperatures. 

Three unknown parameters x,, kp, Ka, in the 
model equations were simulated by minimizing 
errors between experimental data and model 


results. 


4. Modeling results and discussion 


4.1. Estimation of x, 


All values of x,, partition factor of strongly 
adsorbed EYO, are nearly equal to 0.56 by simu- 
lation at 35.0—75.0°C and 25.0—36.0 MPa. It can 
be reasoned that x, is an EYP characteristic, 
which represents the concentration relation be- 
tween weakly adsorbed EYO and strongly ad- 
sorbed EYO. The percentage of strongly adsorbed 
EYO in EYP is) x/Xo x 100% =0.56/0.85 x 
100% = 66%, then the percentage of weakly ad- 
sorbed EYO in EYP is 100 — 66% = 34%. 


4.2. Estimation of K,a, 


In Eqs. (6) and (7), K; is external mass transfer 
coefficient and a, is specific surface area of EYP 
particles. Because a, is difficult to be measured 
accurately, a, and K;, put together were Ky,a,, 
which is called overall volumetric mass transfer 
coefficient. This value is of important for the 
scale-up of this extraction process. The factors 
that influence the value of K,a, are the velocity of 
CO, through the extractor and the size and the 
shape of EYP particle. K,a, at different CO, flow 
rate is shown in Table 3, where experimental 
conditions are 45.0°C temperature, 32.0 MPa 
pressure, 1.2 mm particle size, 250 g of EYP fed. 
In Fig. 5, In(Kja,) is found to have a linear 
relation with CO, flow rate through extractor, Q. 
By linear regression method, the relation between 
K,a, and Q is represented by 


In Ka, = — 4.172 + 0.548 In Q (9) 


Table 3 
Simulation results of Kya, at different flow rates of CO, 





O (kg h~!) K,a, (s~!) RAAD®* 
4.95 0.03776 0.011 
7.19 0.04562 0.012 
7.19 0.04370 0.010 

11.76 0.06030 0.015 


11.91 0.05993 0.015 








* RAAD, Relative absolute averages deviation between ex- 
perimental data and calculated values. 
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Fig. 5. Volumetic mass transfer coefficient (K,a,) vs. flow rate 
of CO, . 


or 
K,a, = 0.015420°4 (10) 


If the flow rate of CO,, Q, is substituted by the 
interstitial velocity of CO, through extractor, u, 
u= Q/S@z<, where S is the area of extractor sec- 
tion, Eq. (10) is replaced by 


K,a, = 1.17u°°** (11) 


The exponent in (10) or (11) lies between the 
value of 0.33 reported for mass transfer liquid 
[20,21] and 0.6 reported for mass transfer in gases 
[22] at comparable Reynolds numbers. Since su- 
percritical fluids have physical and transport 
properties, such as diffusion and viscosity, which 
fall between those of liquid and gases. The expo- 
nent of 0.548 in Eqs. (10) and (11) is plausible. 


4.3. Estimation of adsorption equilibrium 
constants, k,, 


k, is affected by temperature and pressure of 
SC-CO,. Table 4 shows adsorption equilibrium 
constant, k,, simulated by the model. k, increases 
with increase of temperature and decreases with 
increase of pressure. At a constant temperature, 
the density of SC-CO, increases with increase of 
pressure, resulting in solution power increases 
[23,24]. Hence k, decreases with pressure. But it 


Table 4 
Simulated adsorption equilibrium constants k, at different 
temperatures and pressures* 





Temperature (°C) Pressure (MPa) 





25.0 28.0 32.0 36.0 





35.0 76.22 64.52 54.64 48.15 
45.0 76.28 tone: | Sees S31 
55.0 90.09 72.36 57.84 55.93 
65.0 90.91 90.01 81.97 57.14 
75.0 100.00 92.00 73.69 65.62 





* The other extraction conditions are particle size of 1.2 mm 
in diameter, CO, flow rate of about 7 kg h~!. 


seemly violates the principle of decrease of k,, with 
increase of temperature. Following two aspects 
control the solution power of EYO in SC-CO,: 
density of SC-CO, and volatility of EYO, which 
increase resulting in increase of solubility. But at a 
constant pressure, increase of temperature causes 
two controversial effects, i.e. decrease of density 
of SC-CO, and increase of volatility of EYO. As 
a result, solution power of EYO is determined by 
the larger contribution one. In the range of tem- 
peratures used in this work, decrease of SC-CO, 
density as increase of temperature is the larger 
contribution one, causing decrease of solution 
power of EYO in SC-CO,. Therefore increase of 
temperature causes increase of k,. 





4.8 


P=25MPa 
P=28MPa 
P=32MPa 
——_ .  ., 
a | a Y Gis 
1000/T 











Fig. 6. Relation of adsorption equilibrium constants with 
temperature. 
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The adsorption equilibrium constants are also 
shown on a van’t Hoff plot in Fig. 6. The slopes 
obtained by linear regression method at pressures 
of 25.0, 28.0, 32.0, 36.0 MPa are — 758, — 917, 
— 1054, — 753, respectively. The van’t Hoff equa- 
tion is 

AH 

Ink, = In A — =, (12) 
where A is pre-exponential factor, @H is heat of 
adsorption, R is the perfect gas constant, T is the 
temperature. The heats of EYO adsorption on the 
matrix at different pressures are 6.3 kJ mol™! 
(25.0 MPa), 8.2 kJ mol~! (28.0 MPa), 8.8 kJ 
mol~! (32.0 MPa), 6.3 kJ mol! (36.0 MPa). 
Consequently the heat of EYO adsorption is 6—9 
kJ mol~', which belongs to physical adsorption. 
Using a similar method, Goto et al. [15] obtained 
adsorption heat of peppermint oil. 


4.4. Verification of the model 


All parameters of the mathematical model have 
already been determined. Therefore the model can 
be used to predict the concentration of EYO in 
fluid or solid as a function of time or height of 
extraction bed and compare it with experiment for 
verification of the model. Extraction process was 
terminated after a given time. In other wards, the 
CO, flow was then stopped and extractor pressure 
reduced to atmospheric pressure. Aliquot EYP 
samples were abstracted from extractor along 
height of the fixed bed, then EYO concentration 
of each sample was analyzed. The profile of EYO 
concentration of EYP along the bed are shown in 
Figs. 7 and 8, comparing with values predicted by 
the model, which suggests that predicted values 
coincide with the experimental data from the an- 
gle of solid concentration profiles. 

The values calculated by the model at different 
flow rates of CO, are shown in Fig. 9, in which it 
gives good agreement with experimental results. 

Further evidence of the validity of the mathe- 
matical model is provided by Fig. 10, in which the 
feed loaded is changed and recovered oil is shown 
as a function of time. Fig. 10 also represents the 
predicted values by mathematical model in which 
the parameters K,a, (determined by Eq. (10)), x, 





1.0 











1.0 


Fig. 7. EYO concentration of EYP as a function of reduced 
bed length after 38 min. Dots: experimental data; line: pre- 
dicted value. Conditions: temperature, 45.0@; pressure, 32.0 
MPa; CO, flow rate, 12.4 kgh~'; particle size, 1.2 mm in 
diameter; EYP, 0.300 kg. 


and k, are 0.04545 s~', 0.56 and 55.65, respec- 


tively. There is a good agreement between the 
experimental data and calculated values. 


4.5. Application of the model on a pilot plant 


This mathematical model was applied on a pilot 
plant of extraction of EYO with SC-CO,.The 
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Fig. 8. EYO concentration of EYP as a function of reduced 
bed length after 46 min. Dots: experimental data; line: pre- 
dicted value. Conditions: temperature, 45.0@; pressure, 32.0 
MPa; CO, flow rate, 7.2 kgh~!; particle size, 1.2 mm in 
diameter; EYP, 0.300 kg. 
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Fig. 9. EYO recovery versus time obtained by experiment and 
simulation at different flow rates of CO,. Conditions: tempera- 
ture, 45.0@; pressure, 32.0 MPa; particle size, 1.2 mm in 
diameter; EYP, 0.250 kg. 


pilot plant extractor is 5 x 10~* m? in internal 
volume. The pilot plant ran at a temperature of 
45.0°C, a pressure of 32.0 MPa and a CO, flow 
rate of 22.3 kg h~', using a particle size of 1.2 
mm in diameter. The pilot plant was fed 1.650 
kg of EYP. The values of the mathematical 
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Fig. 10. Extraction rate with different amounts of feed. Condi- 


tions: temperature, 45.0@; pressure, 32.0 MPa; particle size, 
1.2 mm in diameter; CO, flow rate, 7.19 kgh—!; 
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Fig. 11. EYO extracted in the pilot plant as a function of time 
and compared with value predicted by the mathematical 
model. Dots: pilot plant data; line: predicted value. 


model parameters K,a, (calculated by Eq. (11)), 
X,, and k, are 0.08453 s~', 0.56, 55.56, respec- 
tively. Both pilot plant experimental data and 
values predicted by the model are all shown in 
Fig. 11, in which recovered oil is plotted as a 
function of time. The model can fairly simulate 
and predict the extraction process in pilot plant. 


5. Conclusion 


According to experiments the controlling steps 
of SC-CO, extraction of EYO from EYP are the 
external mass transfer and the equilibrium be- 
tween two phases. Based on above mechanism, a 
mathematical model was proposed. The parame- 
ters of the model, obtained by simulation, indi- 
cated that: (a) there are 34 wt.% of EYO 
adsorbed weakly on the matrix and 66 wt.% 
strongly; (b) overall volumetric mass transfer co- 
efficient is directly proportional to CO, flow rate 
through extractor with an exponent of 0.548. 
The correlation of the coefficient as the function 
of flow rate of CO, was obtained; (c) the model 
simulated adsorption equilibrium coefficients of 
EYO in the matrix. The heat of adsorption, 
@H, is 6-9 kJ mol™'. 

The mathematical model was verified and ap- 
plied on pilot plant. The results indicated that 
the mathematical model could simulate the ex- 
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traction process of EYO with SC-CO, in a satis- 
factory degree. 


Appendix A. Nomenclature 


specific surface area of EYP, m™' 


particle size of EYP in diameter, m 

EYO recovery, (kg of extracted EYO) (kg 
of EYO in feed)~! 

overall bed length, m 

bed length, m 

overall volumetric mass transfer rate, kg 
m-?s7! 

external mass transfer coefficient, m s~ 
the adsorption constant of EYO between 
the solid and the fluid phases at equi- 
librium, k, = x/y* 

mass flow rate of CO,, kg s~ 
area of extractor section, m? 
time consumed in extraction, s 

superficial velocity of CO, through extrac- 
tor, ms-', U=Q/pS 

fluid interstitial velocity, m s~', u= U/e 
amount of EYP loaded in extractor, kg 
initial EYO concentration of feed, (kg of 
EYO) (kg of EYP)! 

EYO concentration in EYP, (kg of EYO) 
(kg of matrix)! 

strongly adsorbed EYP concentration in 
EYP, (kg of EYO) (kg of matrix)~! 

initial EYO concentration in EYP, (kg of 
EYO) (kg of matrix)~', x) = 0.851 

EYO concentration in SC-CO,, (kg of 
EYO) (kg of CO,)~! 

initial EYO concentration in SC-CO,, (kg 
of EYO) (kg of CO,)~! 

EYO concentration in the solvent in equi- 
librium with EYP (kg of EYO) (kg of 
CO,)~! 

the solubility of EYO in SC-CO,, (kg of 
EYO) (kg of CO,)~! 

reduced bed length, z=h/H 

SC-CO, density, kg m~? 

particle density of the matrix part of EYP, 
kg m~? 


1 


1 


¢ fluid volume fraction 
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Abstract 


The model of broken and intact cells has been applied to the experimental results on oil seed supercritical extraction 
obtained by various authors and on several species of seeds. The results analyzed are related to various extraction 
apparatus ranging from very low laboratory scale to pilot plants. Moreover, a wide range of operating conditions has 
been covered, CO, flow rates ranging between 1.5 and 750 g/min, extraction pressures from 240 to 550 bar, 
temperatures between 25 and 50°C and particle diameters from 0.25 to 4 mm have been considered. Besides, the 
experimental results existing in the literature, systematic scanning electron microscope (SEM) analysis have been 
performed on seed particles belonging to the different seed species studied, thus obtaining data about the microscopic 
cells that bear the oil that are characteristic of the different seed structures. Operating in this manner the differential 
mass balances that characterize this kind of models have been supported by microscopic information on seed 
structure and the number of adjustable parameters in the model has been reduced to only one, the internal mass 
transfer coefficient (4,). A fair good fitting of all the available experimental results has been obtained using best fit k; 
values ranging between 2.4 x 10~’ and 9.2 x 10~* m/s and producing a coherent description of the extraction 
process. Some simulation tests have also been performed that evidenced the role of particle size and of internal mass 
transfer and of their interaction on the overall performance of the extraction process. © 2001 Elsevier Science B.V. 


All rights reserved. 
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fluid processing is very promising for improving 
the traditional industrial techniques that are based 
on hexane extraction for vegetable oils. 

SFE of seed oil has been studied by several 
authors from the processing point of view and a 
wide range of seed species has been explored, 
wheat germ [1]; oats [2]; corn germ [3]; cottonseed 
~ * Corresponding author. Tel.: + 39-089-964116; fax: + 39- [4,5]; soybean [5,6]; rice bran [7]; evening primrose 
089-964057. [8]; jojoba [9]; rapeseed [6]; peanut [10] and grape 

E-mail address: reverch@dica.unisa.it (E. Reverchon). seed [11]. 


1. Introduction 


Vegetable oils are one of the most interesting 
products proposed for supercritical fluid extrac- 
tion (SFE). Indeed, in these cases supercritical 


0896-8446/01/$ - see front matter © 2001 Elsevier Science B.V. All rights reserved. 
PII: S0896-8446(00)00093-0 
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However, despite the relatively large number of 
Species processed, only some models of the SFE 
of seed oil have been published. 

The mathematical models proposed are gener- 
ally based on differential mass balance integra- 
tion. Bulley et al. [12] Lee et al. [13] and Fattori et 
al. [14] assumed that mass transfer resistance was 
only in the solvent phase. In other models, the 
authors considered an internal mass transfer resis- 
tance. King and Catchpole [15] used a shrinking 
core model to describe a variable external mass 
transfer resistance where the solute balance on the 
solid phase determines the thickness of the mass 
transfer layer in the external part of the particles. 

A first attempt at the introduction of a physical 
description of vegetable substrate was proposed 
by Sovova. Indeed, Sovova [16] and Sovova et al. 
[17,18] considered the solid phase as divided be- 
tween broken and intact cells containing the veg- 
etable oil. The merit of this model is to use a 
realistic description of the vegetable structure. Its 
major drawback is the large number of adjustable 
parameters it contains (four in the case of sea 
buckthorn and three in the case of grape seed). At 
least three independent data sets are required to 
proper modeling (for example, three different CO, 
flow rates with all the other process parameters 
set at a fixed value) to obtain from least squares 
analysis the set of parameter values that validates 
the model. Moreover, the authors only tested 
simplified forms of this model. The hypothesis of 
Sovova about vegetable microstructure was ver- 
ified by scanning electron microscope (SEM) for 
almond and fennel particles by our research group 
[19,20]. The internal mass transfer coefficient was 
the only one adjustable parameter remaining in 
the model since all the other parameters were 
calculated from data and from vegetable 
microstructure. 

However, only limited set of data have been 
used until now, that are applicable only to the 
specific material studied and no attempts were 
made to extensively test a mathematical model on 
several seed species. 

Therefore, the scope of this work is to try to 
extend the broken—intact cells model plus mi- 
crostructure analysis to all the seed oils obtained 
by SFE using the data sets available in the litera- 


ture and the SEM image analysis of seed 
microstructures. 


2. General description of the SFE experiments 


The extraction apparatus used for the SFE of 
seed oils usually consists of an extractor (internal 
volume V, internal diameter D,) filled with seed 
particles and one or more separators operated in 
series. The separators are equipped with a device 
that allows the discharge of the extracts at fixed 
time intervals. CO, is delivered by a high pressure 
pump. The instantaneous CO, flow rate and the 
total quantity of CO, used are usually measured 
by calibrated rotameters and test meters, 
respectively. 

When an extraction experiment starts, the ex- 
tractor is pressurized with CO, first, then, when 
the operating pressure has been reached, CO, is 
allowed to flow continuously through the extrac- 
tor at the set flow rate. 

Before extraction, seeds are milled down to 
different mean particle sizes d,. This operation is 
essential for the efficiency of the extraction pro- 
cess — we will see that particle size is one of the 
controlling parameters on the extraction rate and 
on the asymptotic yield of oil. The seed quantity 
loaded in the extractor is indicated by mp and the 
weight flow rate of supercritical CO, is indicated 
by W. The asymptotic yield of oil is the total 
quantity of oil extractable from the seed. 

The main components of vegetable oils are 
triglycerides whose solubility in supercritical CO, 
becomes appreciable at 40°C for pressures around 
200 bar [21]. 

A useful feature that can be added to the 
extraction process is the fractional separation of 
the extracts to isolate seed oil from waxes [22,23]. 


3. Mathematical model 


The modeling of the extraction process is based 
on the following hypotheses. 

1. We suppose that the behavior of all com- 
pounds extracted is similar and can be de- 
scribed by a single pseudo-component with 
respect to the mass transfer phenomena. 
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. Concentration gradients in the fluid phase de- 
velop at larger scales than the particle size. 

. The solvent flow rate, with superficial velocity 
u, is uniformly distributed in all the sections of 
the extractor. 

. The volume fraction of the fluid, ¢, is not 
affected by the reduction of the solid mass 
during extraction. 

Further hypotheses regard the natural 
matrix. 

. The solute in the solid is present in two sepa- 
rate phases. One phase includes the solute 
contained inside the internal structure of the 
particles (the ‘tied solute’ phase). It fills a 
fraction ¢, of the overall volume occupied by 
the seed particles. This value does not change 
during the extraction process and, therefore, 
@, 1S considered constant. The average tied 
solute concentration is called P. The other 
phase is made of the solute freely available on 
the particle surface. The concentration here is 
always the same and, according to our hy- 
potheses, it is equal to the pure solute density 
Po: 

. The fraction of the seed volume filled by the 
free solute before extraction is d= 1 — ,. 

. The fraction of the seed occupied by the free 
solute during the extraction is wd, where is 
Wl. 

. A linear equilibrium relationship applies be- 
tween phases. 

More details on model hypothesis have been 

previously published [19,20]. 

During the extraction, both the free solute and 
the tied solute move from the seed particles to the 
solvent. This flux might reduce the volume occu- 
pied by the seed particles and probably change ¢,; 
and ¢, values. In particular, during the extraction, 
the volume occupied by the free solute phase in 
the untreated particles is filled by the solvent. The 
small value of #, makes negligible the effective 
increase of the volume fraction occupied by fluid 
in the extractor (¢) due to the reduction of the 
volume occupied by the free solute phase during 
the extraction. The inclusion of ¢ and @ changes 
in a mathematical model would make its resolu- 
tion considerably harder; we avoided it according 
to hypothesis d). 


According to the above hypotheses, the mass 
balance on the solute in the extractor is: 
oC 
62 
ow 


—— = () 
Ot wy 
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OP 
+ (1 — Pr" Ps’ a+ (1 — 8) be ° Po’ 


where D, is the axial dispersion, u the superficial 
velocity, p, the fluid density, which we suppose is 
not affected by the presence of the solute, p, the 
bulk density of the non-soluble solid that is the 
mass of non-soluble solids in the seeds per unit of 
filled particle volume, that is the total volume of 
the particle minus the volume of broken cells. 

The general mass balance on the phase of the 
free solute alone is: 


OY  —_—ka(P. — KyC) 


- 0 
Po ar rt es until ¥>0, 





otherwise Co =0 (2) 
where a is the specific surface of the seed particles 
and k, is the external mass transfer resistance. 

We simulated the initial equilibrium regime that 
has been experimentally observed by the authors 
in the literature by using very large values of k; in 
Eq. (2). The correctness of this procedure was 
verified by changing these k, values and observing 
that no significant variations occurred in the mod- 
eled extraction kinetics. 

The mass balance on tied solute is: 


ae oP KO) 
Ot a Ss (1 — &)@, 


where k; is the internal mass transfer resistance. 

In Eqs. (2) and (3) the same mass transfer area 
a was included. This assumption is not strictly 
true since the mass flux might affect different area 
partitions of the particle surface. This model how- 
ever is unable to make a distinction between the 
effect of a mass transfer area change and that of 
a mass flux change. For this reason, the exchange 
area was conventionally put equal to a. 

The system of equations from Eq. (1) to Eq. (3) 
has a unique solution when the initial conditions 
(i.c.) on C, P and w and the boundary condition 
(b.c.) on C are given: 





(3) 
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fic.) At t=0: C=C; P=P,; = P,; 


for each z (4) 


l 
-C-—D,: 
é 


(b.c.) At z=0: 


for each f (5) 


(b.c.) At z=Z: <=0 for each f (6) 

We preliminary performed some tests to evalu- 
ate the influence of the axial dispersion on the 
SFE process. The value of D,; we used was ob- 
tained from the correlations proposed by Fu- 
nazukuri et al. [24] and Catchpole et al. [25] using 
two different experimental methods. These two 
correlations give similar D, values [24]. However, 
the presence of D, practically did not affect the 
modeling results; i.e. the same results were ob- 
tained with and without inserting the axial disper- 
sion. Therefore, we decided to perform modeling 
omitting the axial dispersion in Eq. (1). The corre- 
sponding boundary conditions 4 and 5 were sub- 


stituted by: 
(b.c.) At z= C=0 _ for each ¢ (7) 


We also supposed that the loading time of the 
solvent into the extractor is long enough to enable 
the fluid to reach the equilibrium concentration 
before the extraction started. Therefore, from hy- 
pothesis h) we have: 


Po = KyCo (8) 


where K, is the equilibrium constant between the 
solvent and the free solute phase. This means that, 
before the extraction starts, a part of the oil 
present in the seed saturates the fluid. This condi- 
tion, together with Eq. (8), implies that the tied 
solute phase and the solvent are at equilibrium at 
the extraction start and, therefore, after hypothe- 


sis h): 
Po e- K,Co (9) 


where K, is the equilibrium constant between 
these two phases. As a consequence of this hy- 
pothesis, free solute alone participates to the ini- 
tial saturation of the solvent. This condition is 
necessary for Eq. (8) to apply. Therefore, the 
volume fraction wW, of the free solute at the extrac- 


tion start is different from unit that is the value of 
wy in the untreated seed. 

Both Eqs. (8) and (9) are valid under the hy- 
pothesis of a linear equilibrium relationship be- 
tween the solvent and any of the solute phases. 
This provides a simplified view of the system. 
However, we do not have more accurate informa- 
tion on thermodynamic behavior, to confidently 
assume more complex equilibrium models. 

The value of C, can be calculated from the 
experimental plot of the oil yield as a function of 
the mass of solvent flowed, m,: 


| Mig Mm, Mo 

where m, is the mass of the extracted oil and mz is 
the initial mass of seed. The last term of Eq. (10) 
is the slope of the linear section of the yield curve. 

The values ¢, and ¢, can be determined from 
SEM images of the seed particles and the values 
Y, Po and p, can be determined from the follow- 
ing mass balances. The free oil mass balance that 
relates Cy and Wp to the loading of the untreated 


solid: 


Pl — €)be = p.U1 — €) bho + Copré (11) 
mass balance on the total free oil: 


PeV x = PoPe + Pops, (12) 
mass balance on the non-soluble solids 


ps, = pl — Y.,) (13) 
where Y,, is the value of the asymptotic ex- 
tractable oil yield and p, is the density of seeds. 

From Eqs. (11)—(13) we can calculate ¥%, Po 
and p,, respectively. 

The set of differential Eqs. (1)—(3) was numeri- 
cally integrated using a finite difference method. 
A Wendroff numerical cell was used. This method 
is characterized by an implicit computational cell 
that assures calculation stability whichever is the 
interval (of time and of space) used. We also 
selected the computation intervals to set the 
Courant number equal to | to minimize the com- 
putational error. The solution was performed in 
an explicit manner (calculus simplicity) but the 
method maintains the advantage of the implicit 
cell (calculus stability). This method has been 
described in detail in Lapidus and Pinder [26]. 
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4. Results and discussion 


A scanning of the scientific literature revealed 
that there is enough information (data, experi- 
mental apparatus and vegetable charge character- 
istics) to model SFE for the extraction of 
vegetable oil of sunflower [27], tomato [28], cori- 
ander [29] and grape seed [18]. Almond [19] and 
fennel oil extraction [20] that has been previously 
modeled using the same model have also been 
taken into account. Moreover, liquid SFE extrac- 
tion of peanut seeds [30] was also considered. 
Some other papers that could have been of inter- 
est for our work were discarded since they contain 
not enough information about the experimental 
conditions or about the experimental apparatus or 
the experiments were performed only at different 
pressures, whereas, the other relevant process 
parameters were not varied at all. 

The relevant experimental data for the various 
works in the literature has been summarized in 
Table 1. In the case of tomato oil extraction, the 
value of mp) was not reported in the paper [28]. 
However, we assumed that the volume fraction of 
the fluid phase in the extractor (¢) was 0.4 that is 
a usual value for this parameter and, since the 
extractor volume was given in the paper, we cal- 
culated for m, (1.26 g) as reported in Table 1. 

The second kind of information requested to 
apply the proposed mathematical model is related 
to seed microstructure, particularly to the geomet- 
ric characteristics of the oil-bearing cells. These 
cells should not be confused with biological cells 
that are much smaller. Systematic SEM analysis 
of seed particles has been performed to obtain this 
data. The analysis of seed structures showed that 
the oil-bearing cells of the various seeds studied 
have different shape and size but shape and size 
are fixed for a given species of seed. This last 
observation means that SEM images are required 
only at one time for each seed species. Some 
examples of SEM images of seed particles are 
shown in Figs. 1 and 2 and are related to corian- 
der and grape seed, respectively. In these cases, 
cells are spherical and do not contain any oil since 
images have been produced from SFE processed 
material. However, some starch aggregates are 
still present and form small spherical particles 


located inside the vegetable structure, since this 
last material is not extractable by supercritical 
CO,. SEM images of non-spherical seed cells were 
previously published for fennel seed [20]. 

To evaluate the total volume of the broken cells 
on the particle surface, we first determined their 
number as the ratio between the area of the 
particle surface (spherical particles) and the area 
of the single broken cell, obtained from SEM 
images. From this value, and making simple geo- 
metric considerations on the shape of the cell, we 
obtained the volume of the cell. The total volume 
occupied by broken cells on particle surface is, 
then, approximately one half of the total volume 
of broken cells since only part of these cells is 
available on the particle surface. From these cal- 
culations we obtained ¢, for the various seed 
species and for the used particle sizes, as reported 
in Table 3. 

In Table 2 the characteristic parameters of the 
oil-bearing cells for the various seeds are shown. 

Using the experimental data taken from the 
literature and from microscopic image analysis it 
has been possible to calculate all the parameters 
to be used in the proposed mathematical model, 
except the internal mass transfer coefficient that 
has been used as the only adjustable parameter in 
the model. Table 3 summarizes the values of 
parameters used in the mathematical modeling of 
the various seeds. 

First of all, it is possible to discuss modeling 
results for the single data sets and single seed. 

Sunflower seed oil extraction was performed by 
Perrut et al. [27] using a laboratory and a pilot 
extraction plant. The experiments were performed 
operating at 280 bar and 40°C using CO, flow 
rates ranging from 83.3 to 750 g/min and a parti- 
cle size of 3 mm. The authors modeled the exper- 
imental results assuming that two different 
extraction regimes characterize the process, which 
were related to two differently bonded fractions of 
the solute. They also used a non-linear equi- 
librium curve and found a best fit value for k, of 
2.2 x 10~° m/s. 

SEM analysis of sunflower seed particles 
showed that the oil-bearing cells were ellipsoidal. 
The refit of the experimental data gave us a best 
fit value of k, of 4.0 x 10~7 m/s. The experimental 
data and the model curves are shown in Fig. 3. 
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Coriander seed oil supercritical extraction was 
performed by Catchpole et al. [29]. They did the 
experiments at 250 bar and 40°C using two differ- 
ent CO, flow rates (86 and 213 g/min) and particle 
sizes (0.56 and 0.92 mm). The authors adopted a 
model based on the solid phase mass balance 
performed on a spherical particle (shrinking core) 


Table | 
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using some simplifications that produced an equa- 
tions system with an analytical solution. The best 
fit value of k; they obtained was 9 x 10~°? m/s. 
We used the model with the data reported in 
Tables 1—3 obtaining a best fit value of k, of 
5.4 x 10~*% m/s. The experimental data and the 
model curves were reported in Figs. 4 and 5 where 


Experimental conditions in the literature for the SFE of oil from different seeds 








P (bar)/T (°C) D; (cm) 


H (cm) 


Mo (g) 








Sunflower 
Test 1 
Test 2 
Test 3 
Test 4 


280/40 8.2 28 
280/40 8.2 28 
280/40 8.2 28 
280/40 8.2 28 


Coriander 
Test | 
Test 2 
Test 3 
Test 4 


250/40 10 27 
250/40 10 27 
250/40 10 27 
250/40 10 27 


Grape 
Test | 
Test 2 
Test 3 
Test 4 
Test 5 


280/40 0.8 
280/40 0.8 
280/40 0.8 
280/40 0.8 
280/40 0.8 


Tomato 
Test 1 240/40 0.75 
Test 2 240/40 0.75 
Test 3 240/40 0.75 
Test4 240/40 0.75 
Test5 240/40 0.75 
Test 6 240/40 0.75 


Test 7 240/40 0.75 


Peanut 
Test 1 
Test 2 
Test 3 
Test 4 
Test 5 


Almond 
Test 1 
Test 2 
Test 3 


550/25 
550/25 
550/25 
550/25 
550/25 


300/50 
300/50 
300/50 


Fennel 
Test 1 
Test 2 
Test 3 


200/40 
200/40 
200/40 


500 
500 
500 
500 


1050 
1050 
1050 
1050 


4.5 
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Fig. 1. SEM image of coriander seed particle. Broken cells are 
evidenced. 


the influence of different particle sizes and of 
different CO, flow rates is reported, respectively. 
Very different oil yields were obtained for the 
process times used for the supercritical extraction. 
In the extraction of the larger particles, the most 
inner part of the solid is not efficiently reached by 
the solvent and can participate only in part to the 
extraction process. However, when the extraction 
time is sufficiently long, also in the case of large 
particles the oil yield curves may approach a value 
corresponding to the complete extraction. As we 
will see in the following discussion, this phe- 
nomenon is not limited to coriander seed extrac- 
tion experiments but is a general feature of the 
extraction process when large particles are used. 
Grape seed oil was extracted by Sovova et al. 
[18] operating at 280 bar and 40°C. They studied 


Fig. 2. SEM image of grape seed particle. Broken cells are 
evidenced. 


Table 2 
Shape and characteristic dimensions of the oil-bearing cells for 


the different kind of seeds studied 





Seed Cells shape Mean cell dimensions (1m) 





D, = 60; D, = 10 
D=15 

D=20 

D=25 

L=20; D=18 
D=20 

L=90; D=10 


Sunflower Ellipsoidal 
Coriander Spherical 
Grape Spherical 
Tomato Spherical 
Peanut Cylindrical 
Almond Spherical 
Fennel Hexagonal 





the process behavior for different particle sizes 
(between 0.45 and 0.81 mm) and different CO, 
flow rates (between 0.9 and 1.7 g/min). They 
modeled the experimental results using the broken 
and intact cells model (that, indeed, was first 
proposed by Sovova [16]); however, in their ver- 
sion this model had at least three adjustable 
parameters, since microscopic information was 
not available. The use of five parameters means 
that three or more independent data sets are 
required to obtain a single set of parameter values 
that are the best fit of the experimental data. If 
less than three independent data sets are used 
there will be an indefinite number of parameter 
sets that verify the best fit. However, they ob- 
tained a best fit value of k, of 6 x 10~'° m/s. 

Using the same procedure explained in the pre- 
vious paragraphs, we modeled the experimental 
results on grape oil obtaining a best fit A; value of 
2.4x 10~’ m/s. The model curves together with 
the experimental data are reported in Figs. 6 and 
7 related to the experiments performed at differ- 
ent particle sizes and different CO, flow rates, 
respectively. 

Also in this case the experiments at different 
particle sizes evidenced the yield reduction for 
larger particles in a fixed extraction time (see Fig. 
6) as discussed in the case of coriander seed oil 
extraction. 

Tomato seed oil extraction was studied by Roy 
et al. [28]. They studied the oil yield of the process 
operating at 240 bar and 40°C using various 
particle sizes (between 0.25 and 1.02 mm) and 
different CO, flow rates (between 1.5 and 5.88 
g/min). A high pressure UV—vis spectrophotome- 
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ter was installed at the extractor outlet for the 
continuous monitoring of the extract content in 
the solvent. For this reason, the experimental data 
related to these experiments are practically contin- 
uous and are reported as dotted lines in Fig. 8. 
These authors also proposed the concept of two 
different extraction rates along the process time 


considering broken and unbroken structures. 
They used a model based on differential mass 
balances in which the term related to convective 
mass transfer (in the fluid phase mass balance) 
was substituted with C/t where ¢ is the residence 
time. This simplification transformed the equation 
system from PDEs to ODEs. The best fit value of 


Table 3 
Parameters used in the mathematical modeling of the various seeds studied 














de Co (g/g) Ps (g cm~ *) Po (g/g) 





Sunflower 

Test 1 0.011 0.66 0.42 
Test 2 0.011 0.66 0.42 
Test 3 0.011 0.66 0.42 
Test 4 0.011 0.66 0.42 


Coriander 

Test | 0.004 1.03 0.06 
Test 2 0.004 1.06 0.03 
Test 3 0.004 1.03 0.06 
Test 4 0.004 1.06 0.03 


Grape 

Test 1 0.007 1.00 0.09 
Test 2 0.007 1.00 0.09 
Test 3 0.007 1.05 0.03 
Test 4 0.007 1.04 0.04 
Test 5 0.007 1.01 0.08 


Tomato 

Test 1 0.008 0.81 0.33 
Test 2 0.008 0.81 0.33 
Test 3 0.008 0.81 0.33 
Test 4 0.008 0.99 0.09 
Test 5 0.008 0.92 0.17 
Test 6 0.008 0.86 0.25 
Test 7 0.008 0.81 0.33 


Peanut 

Test 1 0.018 0.64 0.70 
Test 2 0.018 0.68 0.62 
Test 3 0.018 0.78 0.41 
Test 4 0.018 0.82 0.33 
Test 5 0.018 0.91 0.20 


Almond 

Test 1 0.011 0.54 0.98 
Test 2 0.011 0.62 0.73 
Test 3 0.011 0.54 0.98 


Fennel 

Test | 0.002 1.08 0.1 
Test 2 0.002 1.08 0.1 
Test 3 0.002 1.08 0.1 
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Fig. 3. Yield curves of sunflower oil at various flow rate of 
CO,, d,=3 mm; M, W=750 g/min; @, W= 416.6 g/min; A, 
W = 166.6 g/min; V, W= 83.3 g/min. 


k; resulted 4.1 x 10~’ m/s. In this case too the oil 
yield obtained in the used extraction time de- 
creased with the increase of seed particle size (see 
Fig. 9) confirming the conclusions reached in the 
previous paragraphs. 

We modeled tomato extraction yield data in the 
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Fig. 4. Yield curves of coriander oil at various average particle 
diameters, W = 86 g/min; M, d, = 0.56 mm; ®@, d, = 0.92 mm. 
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Fig. 5. Yield curves of coriander oil at various CO, flow rates, 
d, =90.56 mm; @, W=213 g/min; @, W=86 g/min. 


same manner as previously described and ob- 
tained a best fit value of k; of 3.0 x 10~’” m/s. 
Peanut oil was extracted by Goodrum et al. [30] 
using liquid CO, at 550 bar and 25°C. They used 
a flow rate of 72 g/min and particle diameters 
between | and 4 mm. Though liquid and not 
supercritical CO, was used in these experiments, 
we decided to use this data since we expect that 
liquid CO, can show a behavior similar to super- 
critical CO, in this kind of processes. Fitting of 
the experimental results was performed using an 
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Fig. 6. Yield curves of grape oil at various average particle 
diameters, W=1.7 g/min; M, d,=0.45 mm; ®@, d,=0.73 
mm; A, d,=0.81 mm. 
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Fig. 7. Yield curves of grape oil at various CO, flow rates, 
d, = 0.36 mm; M@, W=1.7 g/min; @, W=0.89 g/min. 


empirical equation. Therefore, no value of k; was 
found in this case. 

In this case large particles diameters (tests 2—5) 
give very low values of free oil and very different 
asymptotic oil yields for the given extraction time 
were obtained. This fact confirms the influence of 
particle size on the efficiency of the extraction 
process also when liquid CO, is used. 

Our modeling results have been obtained in the 
usual manner and are reported in Fig. 10 together 
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Fig. 8. Yield curves of tomato oil at various CO, flow rates, 
d,, = 0.25 mm. 
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Fig. 9. Yield curves of tomato oil at various particle diameters, 
W = 3.42 g/min. 


with the experimental data. In this case the best fit 
value of k; is 3.0 x 10~7 m/s. 

Almond oil supercritical extraction was per- 
formed by our research group [19] operating at 
300 bar and 50°C at two different particle diame- 
ters (0.3 and 0.7 mm) and at two different CO, 
flow rates (12 and 23.3 g/min). The same model as 
the one proposed in this work was used with only 
one adjustable parameter. However, we refit the 
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Fig. 10. Yield curves of peanut oil at various average diame- 
ters, W=72 g/min; @, d,=1.1 mm; @, d,=1.4 mm; A, 
d,=2 mm; V, d,=2.8 mm; @, d,=4 mm. 





E. Reverchon, C. Marrone / J. of Supercritical Fluids 19 (2001) 161—175 


experimental results since now we are using a 
more efficient numerical routine. We obtained a 
different best fit value for k; with respect to our 
previous work with k, is 9.2 x 10~* m/s. The 
decrease of oil yield with particle size was ob- 
served also in this case. No figures are proposed 
for almond oil extraction since similar modeling 
results were previously published [19]. 

Fennel oil extraction was also performed by our 
research group [20] operating at 200 bar and 40°C 
using 0.372 mm particle and various CO, flow 
rates (between 8.3 and 24.9 g/min). The model 
used is the one proposed in this work. No refit of 
the experimental data was necessary in this case 
since the same numerical routines were used. The 
best fit value of k; is 8.0 x 10~* m/s. No figures 
are proposed for this case since the same results 
have been previously published [20]. This case has 
been added in the work only for sake of complete- 
ness of the comparison. 

Some general characteristics emerge from an 
overall analysis of the experimental data and of 
the modeling results. 

Solubility (C,) values calculated from the initial 
slope of the experimental data are similar to the 
ones obtained using the Del Valle and Aguilera 
correlation [30]. These authors proposed the fol- 
lowing equation: 


(exp(40.361 — 18708/T + 2186840/T7)*(0.001 p)!°-7?4) 
p 





Ss 


(14) 


to describe the oil solubility in supercritical CO . 
The correlation is valid in a range of pressures 
between 150 and 800 bar. Where C, is the solubil- 
ity of oil in supercritical CO, expressed in g/g, T 
the temperature in Kelvin and p is the density of 
the supercritical fluid in kg/m*. Table 4 shows the 
comparison between the values of Cy and those 
evaluated using Del Valle and Aguilera correla- 
tion [31]. 

As it can be seen from Table 4, the values of Cy 
and C, are in good agreement, the maximum 
deviation between solubility values calculated 
from the experimental data and from correlation 
was obtained for tomato seed oil, where Cy is two 
times larger than C,. We used this correlation also 


Table 4 
Equilibrium oil concentration at different extraction 
conditions* 





Seed Co (g/g) C, (g/g) 





Sunflower 0.011 0.006 
Coriander 0.004 0.005 
Grape 0.006 0.006 
Tomato 0.008 0.004 
Peanut 0.018 0.015 
Almond 0.011 0.007 
Fennel 0.002 0.002 





* Cy has been calculated from the initial slope of the experi- 
mental data, whereas, C, has been evaluated using the Del 
Valle and Aguilera [31] correlation. 


to evaluate peanut oil solubility in liquid CO,. 
Also in this case, the correlation gives a solubility 
value in reasonable agreement with the experi- 
mental data. These results confirm that the seed 
oil located on the particle surface is not connected 
(for example, adsorbed) on the solid structure and 
solubility equilibrium controls the first part of the 
SFE process. These conclusions are in agreement 
with some of the hypothesis of the proposed 
model. 

A second general conclusion is that the use of 
particle diameters larger than about 0.5 mm, in 
connection with internal mass transfer resistance 
values and with the extraction time, produces an 
apparent or real reduction of the asymptotic oil 
yield. This phenomenon has been attributed to a 
low efficiency of the solvent access to the inner 
part of the particles. Indeed, for very long extrac- 
tion times similar asymptotic oil yields should be 
obtained also for very large particle sizes despite 
the low efficiency of the solvent penetration in the 
core of larger particles. When the different sets of 
experiments have been performed at the same 
particle size, as in the case of sunflower (3 mm 
particles), this phenomenon is obviously not 
evidenced. 

Moreover, the experimental results for larger 
particles practically do not show the first linear 
part of the extraction yield. It can mean that the 
quantity of broken cells on the particle surface is 
negligible when compared with the intact cells or 
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the quantity of oil contained in broken cells is not 
sufficient to reach equilibrium conditions. This 
phenomenon is particularly evident in the case of 
peanut oil extraction in which experiments were 
performed using particles with mean diameters up 
to 4 mm (see Fig. 10). 

However, the production of very small particles 
can largely increase the grinding costs and could 
produce bed caking with the formation of chan- 
nels along the extraction bed in which the super- 
critical solvent can preferentially flow, thus 
reducing the extraction efficiency. Therefore, an 
optimization of particles diameter is necessary. 

The limit behaviors of very large particles with 
very reduced quantities of free oil and of very 
small particles with practically all freely available 
oil are correctly described by the proposed mathe- 
matical model. Indeed, if very large particles are 
used the free oil quantity becomes negligible and 
o;=0 and ¥=0. Therefore, Eq. (2) is automati- 
cally excluded from the computation. When the 
seed particles are very smail it results ¢,=0 and 
P=0 and Eg. (3) is automatically excluded from 
the computation. 

We obtained a fair good modeling in all the 
studied cases. Though better fittings have been in 
the case of seeds with the higher oil content. 
However, the comparison of our modeling results 
with those obtained by the referenced authors, 
shows that they also obtained a general fair good 
modeling of the experiments. Indeed, the superi- 
ority of our model when compared with the ones 
proposed by the other authors does not rely on a 
better fitting of experimental results, but its major 
merits descend on the overall coverage of results; 


Table 5 


the consistent physical hypothesis and the consis- 
tent calculation of k; values for the various seeds. 

To support this last consideration in Table 5, 
we summarized the best fit values of the internal 
mass transfer coefficient k; for the various seeds 
obtained in this study and those obtained by the 
other authors. 

The best fit values of k; found for the various 
seeds in this work span along less than two orders 
of magnitude whereas, the values obtained by the 
other authors cover about five orders of magni- 
tude. This result means that a more consistent 
evaluation of k, has been obtained in our study. 

Another point to be discussed is the influence of 
the seed structure on the internal mass transfer 
resistance. The seed structure can be described as 
formed by cellulose and lignin. Lignin is a very 
compact natural polymer with a higher degree of 
polymerization than cellulose. Therefore, lignin 
could be practically inaccessible to the supercriti- 
cal fluid and lignin content of the various kind of 
seed can strongly influence mass transfer resis- 
tance, i.e. the internal mass transfer resistance 
increases with lignin content of the seed. The 
existence of lignin structures can also give another 
explanation for the oil yield reduction for the 
larger seed particles, in these particles some of the 
seed lignin structures could not be destroyed by 
grinding; therefore, there will be some internal 
part of the particle that cannot be accessed by the 
solvent. Consequently, the aliquot of oil con- 
tained into these inaccessible structures is not 
extracted at all and a real oil yield reduction is 
obtained. 


Best fit values of the internal mass transfer coefficient k; for the various seeds proposed in this study and by the other authors that 


performed seed oil supercritical extraction modeling 











Seed 


P (bar)/T (°C) 


k; (m/s) 








Sunflower 
Coriander 
Grape 
Tomato 
Peanut 
Almond 
Fennel 


280/40 
250/40 
280/40 
240/40 
550/25 
300/50 
200/40 


40x 1077 
40x 1078 
2.4x 1077 
3.0x 1077 
3.0x 1077 
92x1078 
8.0x 10~8 








k; (m/s) References 











2.2x 107° 
9x 10~° 
6x 107 !° 
411077 


7.5 x 107? 
8.0x 1078 


Perrut et al. [27] 
Catchpole et al. [29] 
Sovova et al. [18] 
Roy et al. [28] 
Marrone et al. [19] 
Reverchon et al. [20] 
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Fig. 11. Simulation of oil yield curves from a generic seed at 
different values of the particle diameter for k; = 3 x 10~® m/s. 


At this point, we tried to simulate the oil yield 
behavior for a generic seed assuming that super- 
critical extraction has been performed at 280 bar 
and 40°C using a CO, flow rate of 25 g/min and 
that the oil-bearing cells are spherical with a mean 
diameter of 20 um. The oil content is the seed was 
fixed at 40% by weight. Then, we used the mathe- 
matical model for three particle diameters (0.1, 
0.5 and 2 mm) and with two different internal 
mass transfer coefficients, 3 x 10~° and 3 x 10~° 
m/s. The resulting oil yield simulated curves are 
reported in Figs. 11 and 12. 

The first linear part of the oil yield curves is 
showed to depend not only on the free oil extrac- 
tion but also on the extraction of tied oil located 
near the particle surface and is related to the value 
used for the internal mass transfer coefficient. It 
indicates that equilibrium conditions are reached 
also during the extraction of oil from unbroken 
cells located immediately under the particle sur- 
face. Indeed, from the comparison of Figs. 11 and 
12 it is possible to see that very different extrac- 
tion kinetics are obtained when the same particle 
size (and, thus, the same quantity of broken cells) 
is used with different values of k,;. The increase of 
the particle diameter together with the increase of 
mass transfer resistance is also demonstrated to 
control the shape and the attainment of the com- 


plete oil extraction. The increase of particle size in 
Fig. 11 with a k; value of 3 x 10~® produces a 
smooth transition between the two extraction 
regimes and the attainment of a quasi asymptotic 
oil yield for all the considered particle sizes within 
600 min extraction. The use of the same particle 
sizes with a k, value of 3 x 10~°? produce very 
different results. A sharp transition between the 
linear part and the mass transfer controlled part 
of the extraction process can be observed for 0.5 
and 2 mm particles. Also the extraction yields 
obtained in 600 min for the same particle sizes are 
very far from the complete extraction value. 

No inaccessible volume has been introduced in 
the simulations performed. Nevertheless, very dif- 
ferent yields at 600 min have been simulated as 
reported in Fig. 12. Though this result is not 
conclusive about the existence or not of a non-ac- 
cessible volume inside the larger particles, it 
demonstrates that a combination of particle di- 
ameter and mass transfer resistance is sufficient to 
produce lower yields with respect to the effective 
oilseed content for the extraction times selected 
for simulation. 

A final question is, can this model be valid for 
all seed structures existing in the nature? Veg- 
etable structures are not only very complex but 
also plants have been adapted to different condi- 
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Fig. 12. Simulation of oil yield curves from a generic seed at 
different values of the particle diameter for k; = 3 x 10~° m/s. 
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tions. Therefore, it happened that we found an 
exception for seed structure when we studied 
hiprose seed oil supercritical extraction. This seed 
is characterized by a very low oil content and 
SEM analysis of its microstructure showed that it 
did not develop specialized oil-bearing structures. 
In this case, seed oil is contained in very long 
channels that form the vascular structure of the 
seed. The oil is still present in aliquots of free oil 
and difficult to be extracted oil, therefore, a differ- 
ent model was applied to explain all the experi- 
mental evidences [32]. 


5. List of symbols 


specific surface of the solid (1/cm) 

solute concentration in the solvent in 
terms of mass of solute per unit mass of 
solvent (g/g) 


solute concentration in the solvent at the 
beginning of extraction (g/g) 

axial dispersion (m7/s) 

particle diameter (cm) 

mass transfer coefficient from the free oil 
phase to the solvent (m/s) 

mass transfer coefficient from the tied oil 
phase to the solvent (m/s) 

seed oil equilibrium constant between the 
tied oil phase and the solvent (g/g) 

seed oil equilibrium constant between the 
free oil phase and the solvent (g/cm?) 
mass of seed charged (g) 

mass of extracted oil (g) 

mass of solvent used (g) 

solute concentration in the solid in terms 
of mass of solute per unit mass of non- 
soluble solid (g/g) 

solute concentration in the solid at the be- 
ginning of extraction (g/g) 

extraction time (min) 

superficial velocity of the solvent (cm/min) 
extraction yield (—) 

asymptotic extractable oil yield (—) 

axial coordinate in the extractor (cm) 


Greek letters 


é  voidage of the extraction bed (—) 


@, fraction of the particle volume filled by 
the free oil phase (—) 
fraction of the particle volume filled by 
the tied oil phase (—) 
density of the untreated particles (g/cm?) 
solvent density (g/cm?) 
vegetable oil density (g/cm*) 
density of the non-soluble solid (g/cm*) 
ratio of the free oil volume with respect 
to the broken cells (—) 
ratio of the free oil volume with respect 
to the broken cells at the beginning of 
extraction (—) 
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Abstract 


The effects of modifying agents and dyeing conditions of dispersed dyes for PET films and fibers has been 
thoroughly studied. This research reports on the dye incorporation process of non-modified and N,N-dimethylacry- 
lamide modified PET fibers with supercritical CO. Spectral analysis in the infrared region to the modified PET fibers 
showed peaks of the modifier, but the thermal stability of the modified fiber did not show variation when compared 
with the non-modified. The amount of disperse dye incorporated in the PET fibers was determined by UV-Vis 
spectroscopy at 579 nm through the dye extraction using N,N-dimethylformamide. A 2* factorial design was realized 
with the purpose to study the influence of variables in the dye incorporation process as well as their interaction 
effects. The pre-treatment of the PET fibers with N,N-dimethylacrylamide increases the amount of incorporated dye 
3.8 times, in average. The main effect on the dye incorporation was the dyeing time for the modified fibers, but for 
the non-modified fibers, the pressure and the temperature presented analogous main effect © 2001 Elsevier Science 
B.V. All rights reserved. 


Keywords: Poly(ethylene terephthalate); N,N-dimethylacrylamide; Disperse dye; Supercritical CO; Factorial design 








1. Introduction 


Woven fibers of poly(ethylene terephthalate), 
PET, represent the largest percentage among syn- 
thetic fibers in the market [1], and are used as 
pure material or mixed with cellulose. It is utilized 
in the making of different articles, like, curtains, 
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caps and others. In order to have a better appear- 
ance and due to the utilization of these articles, 
addition of dyes in the manufacturing process is 
required. 

Among synthetic fibers, PET has the most com- 
pact and crystalline [2] structure, except for some 
kinds of polyester, of low commercial use. These 
fibers are dyed with disperse dyes, insoluble in 
water (5—30 mg/l), which are applied with dispers- 
ing agents to achieve stable dispersions [2]. 
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Several authors [2—8] have studied the effects of 
modifying agents and dyeing conditions of dis- 
persed dyes for PET films and fibers. According 
to Stinson and Obendorf [6], non-modified PET 
absorbs only 0.4% of water, and does not swell in 
water. This lack of interaction between PET and 
the aqueous dyebath requires either the use of 
high temperature or high-pressure systems. These 
conditions allow the dyeing with disperse dyes to 
achieve industrially acceptable dyeing rates. In 
addition, low molecular weight organic com- 
pounds, called carriers (phenols, amines, aromatic 
hydrocarbons, esters, etc.), which are rapidly ab- 
sorbed, have been used to accelerate the dyeing 
rate [2]. 

Saligram and collaborators [7], report a dyeing 
using sonication in order to form and break small 
bubbles, resulting in a temperature and pressure 
increase at a microscopic level. This cavitation 
Causes agitation, making the dyeing better at low 
temperatures achieving levels comparable with 
dyeing at boiling using carriers. 

Knittel and collaborators [8], showed that syn- 
thetic fibers like PET can be dyed using carbon 
dioxide in supercritical state as solvent, avoiding 
the water pollution and the need of drying with 
excellent firmness levels. 

The use of supercritical fluid in the dye impreg- 
nation of textile materials have been studied due 
to its excellent properties like density, viscosity, 
diffusivity and dielectric constant. The dielectric 
constant of supercritical fluids can be significantly 
changed in the supercritical region and generate 
solvent properties in substances that do not ex- 
hibit such properties in their liquid or gaseous 
state [9]. 

According to Valldeperas [10], dyeing studies 
carried out in polyester and carbonic anhydride at 
140°C and 250 bar, using a series of disperse dyes, 
showed that the dye absorption level and the 
fastness property obtained are comparable with 
the ones of the conventional dyes at high temper- 
atures, without the need of the redox process on 
the batch preparation. Although there are no 
studies concerning the economical viability, its 
ecological implications are evident. Valldeperas 
[10] has shown the advantages of the use of 
supercritical CO, in dyeing processes in compari- 
son with the conventional ones. 


Studies on the solubility of azo and an- 
thraquinone disperse dyes have been published 
and may be used as additional information to the 
studies of PET dyeing in supercritical CO, [11]. 

N,N-dimethylacrylamide has been used in the 
biomaterials area, due to its hydrophilic charac- 
ter, in the obtention of polymeric materials com- 
patible to blood and tissue. The grafting of the 
N,N-dimethylacrylamide monomer in polymeric 
substrates like poly(dimethylsiloxane), to be ap- 
plied as catheter [12], and the obtention of selec- 
tive semipermeable amphiphilic membranes [13], 
used to determine the diffusional characteristics of 
glucose, insulin and albumin, are applications in- 
dicating that the N,N-dimethylacrylamide is not 
harmful to human body. 

In previous work [14,15] we showed that the 
pre-treatment of the PET films with N,N- 
dimethylacrylamide, favors the dye incorporation 
in the fiber during the dyeing in an aqueous 
medium. The use of this pre-treatment decreases 
the dyeing time and increases the dyeing incorpo- 
ration in approximately three times, in compari- 
son with the non-modified PET dyeing process. 

The objective of this study was the investigation 
of the supercritical CO,-aided dye incorporation 
process in non-modified and N,N-dimethylacry- 
lamide modified PET fibers. By using a fatorial 
design, the variables of time, temperature and 
pressure of dyeing, as well as treatment with 
modifier, were investigated. 


2. Experimental 


2.1. Materials 


The commercial dye used was Samaron HGS 
(Dy Star)® (Disperse Blue 79) (C,3;H,;BrN,0,). 
The solvent (modifier) was N,N- 
dimethylacrylamide(C;H,NO) (Aldrich). Com- 
mercial samples of PET fibers were kindly 
supplied by COCAMAR, Maringa, Brazil. N,N- 
dimethylformamide was also used to extract the 
dye from the fiber after dyeing process. Fig. 1 
shows the polymer, dye and modifier structures. 
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Fig. 1. (A) Dye, Samaron HGS (Dy Star)" (C I Disperse Blue 


79); (B) Modifier, N,N-dimethylacrylamide; (C) Polymer, 
Poly(ethylene terephtalate). 


2.2. Procedures 


PET fibers in coil form were immersed in N,N- 
dimethylacrylamide modifier solvent, under agita- 




















tion at 90°C for 15 min. Afterwards the modified 
fiber was centrifuged to eliminate the excess of 
N,N-dimethylacrylamide. In the dyeing experi- 
ments, it was used 3 wt.% of dye in relation to the 
mass fiber. The dyestuff amount is then in large 
excess compared with average solubility of dis- 
perse dyes in supercritical CO, [16—18]. The ex- 
cess of dyestuff amount is in the range 4—95 times 
compared with the total amount absorbed from 
the PET fiber. 

Modified and non-modified PET fibers were 
placed in a dyeing cell with supercritical CO, as 
schematized in Fig. 2. CO, coming from the stor- 
age tank (A) is compressed by a piston (B) until 
reaches the requested pressure. The gas goes 
through: a solenoid valve (C); a one-way valve 
(D); a flow distributor (E) and reaches the infu- 
sion cell (F) in the thermostatized bath; (G) main- 
tained at the requested temperature. By using two 
thermocouples (one in the bath and the other in 
the infusion cell) the requested temperature in the 
infusion cell inner region is controlled. There is a 


























Fig. 2. Supercritical fluid equipment scheme; storage tank (A); piston (B); solenoid valve (C); one way valve (D); flow distributor 
(E); infusion cell (F); thermostatized bath (G); pressure valve (H) and manometers (I). TM, temperature measurer; TS, temperature 


set. 
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Table | 


Codified variables of each experiment of the 2* factorial design 





Test Dyeing temperature* Dyeing pressure” 


Dyeing time‘ Pre-treatment with DMACr? 
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* Dyeing temperature (°C), (50 (—); 65(+). 

> Dyeing pressure (kgf cm~ 7), (140(—); 250(+)). 
© Dyeing time, (30 min.(—); 3 h(+)). 

4 Pre-treatment, (without(—); with (+)). 


pressure reducer valve in the exit (H). A manome- 
ter (I) measures the pressure drop in the line 
before the one-way valve. The sample and the dye 
were placed in individual supports inside the infu- 
sion cell. After, the temperature of bath was 
adjusted to the requested value and the exit valve 
was shut. With the compressor already regulated, 
the entrance valve is opened and the timer is 
trigged after the infusion cell has reached the 
requested pressure and temperature. After this 
operation, the exit valve is opened and the CO, is 
liberated to the atmosphere. 

A factorial design was used with the purpose of 
analyzing the influence of each variable studied in 
the dye incorporation process as well as their 
interaction effects. The variables studied are the 
following — dyeing time, dyeing temperature, 
and pressure of CO, — in the dyeing process. 
Each variable was maintained in two levels. Table 
1 shows the codified variables and the levels of the 
factorial design. 


2.3. Measurements 


Differential scanning calorimetry (DSC) mea- 
surements were performed in order to determine 
the glass transition temperature of modified and 
non-modified PET fibers. Thermogravimetric 
analysis (TGA) of non-modified and WN,N- 
dimethylacrylamide modified PET fibers were per- 
formed in order to verify the thermal stability. 
The DSC and TGA measurements were made 
using a Shimadzu equipment, model 50. IR-ATR 
spectra of non-modified and N,N-dimethylacry- 
lamide modified samples were obtained in the 
infrared region using a FTIR Bomen equipment. 
Scanning electron microphotographs (SEM) were 
taken using a JEOL JSM-6340F field emission 
scanning electron microscopy (FESEM). 

The amount of absorbed dye on non-modified 
and N,N-dimethylacrylamide modified PET fibers 
was determined by UV-Vis spectroscopy at 579 
nm{15], through the dye extraction using N,N- 





W.L.F. Santos et al. /J. of Supercritical Fluids 19 (2001) 177-185 18] 


dimethylformamide (DMF) and posterior mea- 
surement of dye absorption intensity. The samples 
were weighted, and then immersed in an extrac- 
tion cell with N,N-dimethylformamide at 80°C. 
The extraction was carried out in three steps of 30 
min each, using, respectively, the following sol- 
vent volumes 20, 20 and 10 mil. 


3. Results and discussion 


Fig. 3 shows FTIR spectra of PET fibers before 
(A) and after immersion in N,N-dimethylacry- 
lamide (B). Spectral analysis of modified PET 
fibers, showed peaks of the modifier at 1647, 1612 
and 1175 cm~' (Fig. 3B) which were, respectively, 
attributed to carbonyl stretching, C=C stretching 
and C-—N stretching of the modifier [19]. By com- 
parison of spectra (A) and (B), the absorption 
bands can be attributed to N,N-dimethylacry- 
lamide showing an evidence of it incorporation in 
the PET fibers. These peaks indicate the absence 
of a chemical interaction between the modifier 
and PET, as already observed by Ershov and 
collaborators [20] in studies using polyester fibers 
and acrylic copolymers. 


The analysis of DSC data from non-modified 
(Fig. 4A) and N,N-dimethylacrylamide modified 
PET fibers (Fig. 4B) show that melting tempera- 
ture (253°C) is not changed with the treatment 
process. A better defined glass transition tempera- 
ture at about 70°C was observed in the curve of 
the modified PET fibers (B). The same behavior 
was found by Rubira and collaborators with PET 
films [15]. 

The decomposition temperature, at a 10% level, 
determined by TGA, for non-modified (Fig. 5B) 
and modified PET fibers (Fig. 5B) was 387 and 
379°C, respectively. It is an evidence that the 
treatment does not alter the thermal stability of 
the polymer. 

The scanning electronic microphotographs of 
the non-modified (Fig. 6A) and N,N-dimethy- 
lacrylamide modified PET fibers (Fig. 6B) show 
alterations in the modified fiber surface due to the 
N,N-dimethylacrylamide treatment. The modified 
fiber shows flatter surface than the non-modified 
fiber, probably due to the surface corrosion, but it 
was not possible to quantify the corrosion level. 

The amount of dye absorbed is shown in Table 
2. In order to highlight the data, they were di- 
vided in two groups. One of the groups is related 
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Fig. 3. FTIR spectra of non-modified (A) and N,N-dimethylacrylamide modified (B) PET fiber (obtained in KBr pellets). 
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Fig. 4. DSC curves of non-modified (A) and N,N-dimethylacrylamide modified (B) PET fiber. 
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Fig. 5. TGA curves of non-modified (A) and N,N-dimethylacrylamide modified (B) PET fiber. 


to the factorial design 2° of the N,N-dimethy- 
lacrylamide modified fiber that corresponds to the 
experiments designed as test 1—8 in Tables 1 and 
2. Another group is related to eight experiments 
carried out with non-modified PET fiber, de- 
scribed as test 9-16 in Tables 1 and 2. 


Table 3 shows the main effect and the interac- 
tion effect of the variables involved in the facto- 
rial design, described as test 1-8, for the 
N,N-dimethylacrylamide modified fiber in Tables 
1 and 2. From data in Table 3 can be observed 
that the dyeing time presents a higher main effect 
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on the amount of incorporated dye. A significant 
interaction effect between dyeing time and temper- 
ature variables is also observed. These data match 
very well with the dyeing conditions data 
from literature [4,5,21—23] and can be understood 
considering that when the N,N-dimethylacrylamide 
is used to modify PET fiber, the surface as well as 
the fiber bulk were in perfect dyeing condition as 
pointed out by Needles and Walker [21], Bendak 
and El-Marsafi [4,5] and Sicardi and collaborators 
[23]. For example they showed that even though the 
surface of PET was very crystalline, the interior 
regions are less crystalline [18] and they also 
showed that the treatment of the fiber with organic 
solvents such as dimethylformamide decreases the 
crystalline appearance of the surface, making the 


Fig. 6. Scanning electron microphotographs of non-modified 
(A) and N,N-dimethylacrylamide modified (B) PET fiber. 
Magnification, 3000 x . 


dyeing easier [3—5]. Sicardi and collaborators [23] 
showed that the diffusion coefficients and solubility 
of disperse dyes in PET, during impregnation with 
a supercritical fluid, increase with increasing tem- 
perature. They attributed this result to the increase 
in flexibility of the polymeric chains due to the 
temperature. They also found that rising the pres- 
sure results in a small increase in the diffusion 
coefficient and these findings were justified based 
on the CO, plasticizing effect. Based on our results 
and from data of the literature [4,21,23] it can be 
emphasized that the modifier N,N-dimethylacry- 
lamide plasticize and swell the PET fibers increas- 
ing the chain flexilbility promoting the SCF-ass- 
isted dye infusion. In this way, it should be empha- 
sized that the amount of dye that can be incorpo- 
rated in a PET fiber is limited by the dyeing time 
and longer dyeing times increase the dye uptake of 
the fibers. 

On the other hand, to have a good dyeing of a 
non modified PET fiber we need, at first, to achieve 
the optimum SF-CO, dyeing conditions. Table 4 
shows the main effect and the interaction effect of 
the variables involved in the factorial design, de- 
scribed as test 9-16, for the non-modified fiber in 
Tables 1 and 2. The results show a main analogous 
effect of the pressure and temperature on the 
amount of incorporated dye. These results are 
consistent with the data from other systems in 
which a semicrystalline polymer is annealed by 
plasticization with CO, allowing polymer chains to 
reorient to more thermodynamically favorable po- 
sitions, forming crystallites [24—27]. For example, 
Lambert and Paulaitis [26] and Baldwin and Suh 
[27] observed maxima in the time-dependent sorp- 
tion studies of CO, in PET. The decrease in CO, 
solubility and absorption at long times is attributed 
to an increase in crystallinity of PET. 

By comparing the data of Table 2, it was shown 
out that when the fibers are submitted to the 
N,N-dimethylacrylamide treatment (test 1-8), it 
incorporates, On average, 3.8 times more dye than 
the non-modified PET fibers (test 9-16). From the 
above discussion it may pointed out that the 
treatment of PET fiber with N,N-dimethylacry- 
lamide decreases the crystallinity of PET. The 
dyeing of such fibers using CO, in supercritical 
conditions is easier than the non modified PET 


fibers. 
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Table 2 


Ratio of mass of absorbed dye and mass of fiber according to the 2* factorial design described on Table | 





Test Mass of absorbed dye (M.,) (g) x 10* 


Mass of fiber (M,) (g) Ratio M,/M, x 10% 





0.649 
YB 5 | 
2.243 
2.635 
0.466 
3.001 
Pe SF - 
11.338 
0.309 
0.518 
0.309 
0.701 
0.832 
0.597 
2.008 
2.583 
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0.0983 
0.1122 
0.1072 
0.1164 
0.0803 
0.0989 
0.1131 
0.1343 
0.1059 
0.1150 
0.0990 
0.1208 
0.2133 
0.1192 
0.1329 
0.1228 





Table 3 
Main effect and interaction effects of the variables involved on 


the factorial design, described as test 1-8, in Tables | and 2 





Variable Main effect 





Dyeing time 2zTx i? 

Dyeing pressure 2.1 x 10~+ 

Dyeing 1.7x 10-3 
temperature 


Interaction effects 

Dyeing Dyeing Dyeing 
time x dyeing _ time x dyeing pressure x dyeing 
pressure temperature temperature 

0.70 x 10-3 1.8x 10-3 1.3x 107 





4. Conclusion 


The pre-treatment of the PET fibers using N,N- 
dimethylacrylamide as modifier increases the 
amount of incorporated dye 3.8 times, in average. 
Such process does not present damages to the 
environment, once (i) the modifier is a non toxic 
agent and it can be easily recovered; (11) the 
dyeing process does not produce liquid effluents, 
but liberates only CO, as waste product. In this 
way, the dyeing of PET fibers, using supercritical 
CO,, is an alternative technique for the incorpora- 
tion of dispersed dye in PET. 


Table 4 

Main effect and interaction effects of the variables involved on 
the factorial design, described as test 9 to 16, in Tables 1 and 
9) 


<= 





Variable Main effect 





Dyeing time 0.28 x 10-3 
Dyeing pressure 0.71 x 1073 
Dyeing 0.72 x 107-3 
temperature 
Interaction effects 
Dyeing Dyeing Dyeing 
time x dyeing _ time x dyeing pressure x dyeing 


pressure temperature temperature 
0.15x 10~3 0.70 x 10~3 0.64 x 10-3 
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Abstract 


Solubilities of carbon dioxide in poly(vinyl acetate) (PVAc) were measured at temperatures from 313.15 to 373.15 
K and pressures up to 17.5 MPa. Diffusion coefficients of carbon dioxide in PVAc were also measured at 313.15 K 
and pressures up to 7 MPa. Solubilities and diffusion coefficients of carbon dioxide in molten polystyrene (PS) were 
studied at temperatures from 373.15 to 473.15 K and pressures up to 20 MPa. An apparatus using a magnetic 
suspension balance (MSB) was constructed for the measurements. The solubilities in the PVAc and the PS were in 
good agreement with literature data. The solubility in both polymers were correlated with the Sanchez and Lacombe 
equation of state to within an average relative deviation of 3.6 and 1.6% for PVAc and PS systems, respectively. The 
diffusion coefficients in PS were correlated with free volume theory of Kulkarni and Stern to within 10% of relative 
average deviation. © 2001 Elsevier Science B.V. All rights reserved. 


Keywords: Solubility; Diffusion coefficient; CO,; Polystyrene; Poly(vinyl acetate) 





1. Introduction 


Polymeric foams are widely used as heat insula- 
tors, food trays, and support materials, because of 
their many advantageous characteristics such as 
low thermal conductivity, light weight, and high 
impact strength. CFC-11 (trichlorofluoromethane) 
and CFC-12 (dichlorodifluoromethane) had been 
used as physical blowing agents for these types of 
industrial foaming processes. However, since pro- 
duction of these gases has been prohibited, new 
environmentally benign alternative blowing agents 
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such as carbon dioxide and nitrogen are being 
considered. For effective process design, solubility 
and diffusivity of these gases in polymers are 
necessary. 

Solubility and diffusivity of carbon dioxide in 
molten or thermally softened polymer systems 
have been investigated by several authors. Newitt 
and Weale [1] studied the solution and diffusion 
of hydrogen, nitrogen, carbon dioxide, and ethyl- 
ene in polystyrene (PS) at temperatures up to 463 
K and over the pressure range of 8-30 MPa. 
Durrill and Griskey [2,3] investigated a number of 
gas + polymer systems at temperature of 461.5 K 
Sato et al. [4,5] studied the solubility of carbon 
dioxide and nitrogen in PS, polypropylene, high- 
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density polyethylene at 313-473 K and up to 20 
MPa. Takishima et al. [6] investigated the solubil- 
ity and swelling caused by dissolution of CO, in 
poly(vinyl acetate) (PVAc) at 313 and 323 K and 
up to 9 MPa. 

The pressure decay method is popular, because 
the apparatus is simple and construction cost is 
not expensive. However, it is difficult to apply at 
high temperatures especially for polymer melts, 
because pressure sensors, which have suitable ac- 
curacy and a small inner volume, are not avail- 
able. Furthermore, in the pressure decay method, 
a large amount of the sample (about 5 g) is 
needed, which results in long measurement times 
in the molten state. The electrobalance method, 
which has many advantages such as a small 
amount of sample (short measurement times) and 
a high sensitivity, is also popular in the solubility 
measurement at near room temperature [7,8]. In 
the electrobalance method, several researchers 
used a microbalance placed in a pressure vessel. 
For this case, however, measurements were lim- 
ited to temperatures below 125°C due to the 
microbalance operating conditions. New type 
gravimetric apparatus, in which temperatures of 
the balance and sample can be controlled inde- 
pendently, have been proposed for solubility mea- 
surements at high temperatures [9,10]. These type 
apparatus, however, can only be used for gases 
whose density is such that natural convection 
effects are tolerable. Study of polymer melts with 
dense carbon dioxide in its supercritical state 
(op. =467 kg/m*) would be expected to present 
experimental difficulties. Kleinrahm and Wagner 
[11] developed a unique balance, so-called a mag- 
netic suspension balance (MSB), for accurate 
measurements of fluid densities. The MSB has 
many advantages similar to microbalances except 
that the sample and the balance are isolated. This 
type instrument is expected to be suitable for 
measurement of the gas solubility and diffusivity 
in polymer at high temperatures and high-pres- 
sures. Loockemann et al. [12] and Schnitzler and 
Eggers [13,14] developed the apparatus using the 
MSB for the solubility and diffusivity of gas in 
solid polymer. In this work, a new apparatus 
based on the MSB was developed for measure- 
ments of solubility and diffusion coefficient of gas 


in molten polymer and solubilities and diffusion 
coefficients of carbon dioxide in PS and PVAc 
were measured to confirm the applicability. Fur- 
thermore, new data of diffusivity of carbon diox- 
ide in PS were studied. 


2. Experimental 
2.1. Materials 


Polystyrene (PS, T, = 381.4 K, Mw = 3.30 x 
10°, Mn=1.07 x 10°) was supplied by Asahi 
Chemicals Inc. (Kawasaki). Poly(vinyl acetate) 
(PVAc, T,=303 K, Mn=1.67 x 10°) was pur- 
chased from Aldrich. All the characteristics of the 
polymers were reported by the suppliers. Carbon 
dioxide (>99.5% purity) was obtained from 
Iwatani Industrial Gases Corp. (Hiroshima). All 
chemicals were used as received. 


2.2. Apparatus and methods 


Solubility and diffusion coefficient of carbon 
dioxide in the polymers were measured by a MSB. 
Schematic diagram of the MSB [11] (Rubotherm 
Prazisionsmeftechnik GmbH, Germany) is shown 
in Fig. 1. An electronically controlled magnetic 
suspension coupling is used to transmit the mea- 
sured force from the sample enclosed in a pressure 
vessel to a microbalance. The suspension magnet, 
which is used for the transmitting the force, con- 
sists of a permanent magnet, a sensor core and a 
device for decoupling the measuring-load. An 
electromagnet, which is attached at the underfloor 
weighing hook of a balance, maintains the freely 
suspended state of the suspension magnet via an 
electronic control unit. Using this magnetic sus- 
pension coupling, the measuring force is transmit- 
ted contactlessly from the measuring chamber to 
the microbalance, which is located outside the 
chamber under ambient atmospheric conditions. 

The distinct merit of the MSB is that the mi- 
crobalance can be tared and calibrated during 
measurements. This zero-point correction and cal- 
ibration of sensitivity of the microbalance is im- 
portant for measuring the gas _ solubilities 
accurately, because long sorption times are gener- 
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ally needed. The MSB offers the possibility of 
lowering the suspension magnet in a controlled 
way to a second stationary position (Fig. lb) a 
few millimeters below the measuring position. At 
this position, a sample basket is set down on a 
support by the suspension magnet’s descent, and 
the sample is decoupled from the balance. This 
so-called ‘zero-point position’ allows a taring and 
calibration of the balance at any time, even while 
recording measurements. 

In this work, the MSB was used for measure- 
ment of gas solubility and diffusivity in polymer. 
The MSB can be used at pressures up to 35 MPa 
and temperatures up to 523 K. Resolution and 
accuracy of the microbalance (Mettler AT26l, 
Switzerland) are 0.01 mg and 0.002%, respec- 
tively. The amount of dissolved gas, W, (g) was 
obtained following equation considering gas 
buoyancy: 


W, = W,(P, T) — W,(0, T) 


+ p(P, TIVAP, T)t1 + SP, T)} + Vo] 
(1) 


where W; and W, are readout of the balance at 
temperature, 7 and pressure, P and at zero pres- 
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sure, respectively, p is the gas density (g/cm*) and 
was obtained from the equation of state of Angus 
et al. [15]. The Vp and the V, (cm*) are volumes 
of the polymer sample and a sample basket (in- 
cluding a Petri dish, a sample holder, and the 
measuring load decoupling device, about | cm*), 
respectively. V, was determined by buoyancy 
measurements with an accuracy of + 0.06% using 
high-pressure CO, and N, at temperatures of 
373.15, 423.15, and 473.15 K and pressures up to 
20 MPa. The Vp was obtained from sample mass 
and specific volume, vp (cm*/g). The vp was ob- 
tained from the Tait [16] parameters reported by 
Rodgers [17] for PVAc. The Tait parameters for 
PS were determined in this work using literature 
data of Zoller and Walsh [18] as follows: 


P 
Up = vl 1 — 0.0894 n| 1+ alt (2) 


where 
vy = 0.7884 exp(5.790 x 10-47) (3) 
B(T) = 887.2 exp( — 4.323 x 10~3T) (4) 


The degree of swelling of polymer, S,, is caused 
by gas sorption. Since it is difficult to measure S,, 
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Fig. 2. Schematic diagram of solubility and diffusivity apparatus. 


of molten polymer, we used predicted value of S,, 
by the Sanchez—Lacombe equation of state (S—L 
EOS) [19,20]. S,, is defined as follows 


(1 + S)v(P, T, S) 
v(P, T, 0) 


where S (g-gas/g-polym) is the solubility and v 
(cm3/g) is the specific volume of the polymer 
calculated by the S—L EOS. Since S,, is a function 
of solubility, it is necessary to solve Eqs. (1) and 
(5) simultaneously. Sato et al. [21] reported that 
the S,, predicted by the S—L EOS agreed with the 
experimental data for CO, + PVAc [6] to within 
an average relative deviation of 3.5%. 

Mutual diffusion coefficients of gas in polymers 
were determined by measuring W, versus time 
curves during gas dissolution. We used following 
four assumptions to determine the diffusion co- 
efficient. First, we assumed Fick’s second law. 
Second, the diffusion coefficient D was treated as 
being independent of gas concentration during 
each stepwise gas dissolution. Third, the diffusion 
was assumed to be one-dimensional in the poly- 
mer sample sheet with all diffusing substances 
entering through a single face plane. Fourth, the 
thickness of the polymer L was assumed to be 
constant during the sorption. The appropriate 
solution of the diffusion equation may be written 
[22] as 


(5) 





Sy, = 


W(t) — WO) 





W,(cc) — W,(0) — 


—(2n +1)? 


41? ©) 





8 :07| am] 
m2, exp 

where W, (t) is amount of gas in polymer at time 
t. Polymer sample was mounted in a Petri dish (¢ 
25 x 4 mm) to keep the sample form constant by 
melt molding the sample in the Petri dish. It is 
important to keep the sample form constant dur- 
ing sorption for the diffusivity measurements. 
This leads to reliable mathematical analysis of the 
results. The diffusion coefficient was determined 
from fitting sorption curves during about 0.5 MPa 
stepwise pressure changes to minimize the effect 
of concentration dependence on diffusion coeffi- 
cient. L was obtained by dividing the sample 
volume corrected with S, at average value be- 
tween initial and final pressures by area of the 
Petri dish. 

The schematic diagram of the apparatus is 
shown in Fig. 2. The apparatus consisted of a 
section for high-pressure gas supply, a section for 
mass measurement, and a data acquisition system. 
Pressure of gas was generated by means of a 
high-performance liquid chromatography (HPLC) 
pump (Jasco 880PU) for condensable gas and a 
pressure generator (High Pressure Equipment 87- 


Loa 
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6-5) for uncondensable gas. The pressure was 
measured with Paroscientific pressure transducer 
(46KR, 41.4 MPa F.S., accuracy 0.01% F-.S.). 
Tubes, the high-pressure generator, and the pres- 
sure transducer were thermostated to avoid con- 
densation of the gas. The sorption cell, in which 
the sample basket was enclosed, was surrounded 
by an aluminum block. The aluminum block was 
mounted with a band heater and was surrounded 
with heat insulation. Temperature of the sorption 
cell was measured with an accuracy of +0.05 K 
using a Pt 100 Q RTD probe inserted in the 
sorption cell wall as well as a thermometer (Hart 
Scientific 1502) and was controlled with a PID 
controller and the heater. The RTD probe as well 
as the thermometer was calibrated against an 
intelligent standard thermometer (Kaye Instru- 
ments, X0860, accuracy, + 0.01 K). 

A sample mounted in the apparatus was de- 
gassed about 10 h in vacuum at a given tempera- 
ture, then carbon dioxide was introduced into the 
sorption cell to get the solubility or diffusivity. 
The readings of the balance, temperature, and 
pressure were acquired by a computer until equili- 
bration was accomplished (it took about 10 min 
to 4 h depending on the diffusivity). Subsequent 
measurements were repeated to get isothermal 
solubility data. Pressure change of carbon dioxide 
was carried out stepwise. The pressure changes 
were about 2—4 MPa for the solubility measure- 
ments and about 0.5 MPa for simultaneous mea- 
surements of the diffusivity and solubility. 

The mass of polymer used in the experiments 
was 0.1 to 0.2 g for PVAc and about 0.35 g for 
PS. Stability of the MSB depends on density 
fluctuation, because fluid buoyancy affects read- 
out of the MSB. Fluctuations of the mass, tem- 
perature, and pressure were less than +0.13 mg, 
+0.025 K, and +0.003 MPa, respectively. 


3. Results and discussion 


3.1. Solubility 


The solubility of CO, in PVAc is shown in Fig. 
3 and Table 1. The solubilities increased with 
increasing pressure and decreased with an increase 


in temperature. These temperature and pressure 
dependence of the solubility were usually ob- 
served in condensable gas and polymer systems. 
The solubilities obtained in this work were com- 
pared with literature data. Our results were in 
good agreement with experimental data obtained 
by using an electrobalance method [23] to within 
an average relative deviation of 4.0% and to 
within 1.6% for measurements [24] made using a 
volumetric method [21]. The solubility of Tak- 
ishima et al. [6] at 313.2 K agreed with our results 
to within 8.5% average relative deviation. 

Error analysis was carried out for CO, + PVAc 
system. Uncertainty of solubility was estimated 
less than 4% at 313 K and at pressures between 
0.5 and 7 MPa. However, this uncertainty does 
not include S,, error, because it is hard to estimate 
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Fig. 3. Solubility of carbon dioxide in PVAc, (a) Takishima et 
al. (1990), W, 313.2 K; (b) Matsuura (1998), W, 313.2 K, 
Fukusako (1997), @, 313.2 K, A, 333.2 K, M@, 353.2 K. 
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Table 1 
Experimental solubility of carbon dioxide in PVAc 








Temperature (K) Pressure (MPa) Solubility without S,, correction 


(g-gas/kg-polym) 


Swelling S,, 
(S—-L EOS) 


Solubility with S,, correction 
(g-gas/kg-polym) 








313.15 5.49 
14.55 
27.80 
41.84 
65.92 
81.63 
98.53 
115.99 
129.44 
145.78 
166.11 
189.95 
209.92 
234.51 
257.33 
283.40 

53.93 
100.34 
150.18 
195.87 
241.96 

38.47 

73.70 
107.63 
139.42 
166.01 

29.46 

86.55 
134.41 
162.58 
167.00 


5.51 0.006 
14.65 0.016 
28.20 0.031 
42.75 0.047 
68.22 0.075 
85.22 0.094 

103.87 0.115 
123.54 0.137 
139.01 0.154 
158.19 0.176 


182.74 0.204 
212.59 0.238 


238.65 0.268 
aikeae 0.306 


305.68 0.345 
346.92 0.393 


56.08 0.064 
108.55 0.124 
171.37 0.197 
239.26 0.277 
332.41 0.387 

39.92 0.047 

79.59 0.094 
121.79 0.144 
167.24 0.199 
214.30 0.256 

30.51 0.037 

97.89 0.120 
170.00 0.209 
237.68 0.293 


299.25 0.369 














the uncertainty of swelling, S,, predicted by S—L 
EOS. Hence, we assumed that S,, had a relative 
uncertainty of 5%. Uncertainty of the solubility 
including the S,, error was estimated to be 6.0% at 
7 MPa and 313.15 K. At 17 MPa and 373.15 K, 
the uncertainties for neglecting and including 
swelling, S,,, error were estimated to be a maxi- 
mum of 7 and 20%, respectively. The S,, value 
affects solubility determination considerably. 

The solubility of CO, in PS is shown in Fig. 4 
and Table 2. Although PS was in its glassy state 
below its glass transition temperature of 381.4 K 
at atmospheric pressure, dual mode sorption be- 
havior [7] was not observed at 373.15 K and over 
2.068 MPa. In this experimental temperature and 


pressure range, PS was in its rubbery state, be- 
cause the glass transition concentration of CO, at 
373.15 K estimated with Chow’s method [25] was 
4.84 g-gas/kg-polym. In the Fig. 4, data of Sato et 
al. [4] obtained by using a pressure decay method 
at 373.2 K is also shown and agrees well with 
those of this work. 

The swelling, S,,, predicted by the S—L EOS are 
given in Tables 1 and 2. The S,, showed a linear 
dependence on solubility for both the systems. 
High solubilites result in large degrees of swelling 
in the PVAc system. Although the swelling is not 
so high in the PS system, the relative S,, correction 
of the PS system was comparable with the PVAc 
system at the same temperature and pressure. 
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Fig. 5 shows the temperature dependence Hen- 
ry’s constant, Kp (kg-MPa/cm? (STP)) for CO, in 
PS and PVAc obtained in this work along with 
literature data. The reciprocal Henry’s constants, 
1/Kp in PS of this work agreed with previous data. 
Data for PVAc of Takishima et al. [6] were 4.6% 
higher than those of this work. Since literature 
data of Matsuura [23] and Fukusako [24] were 
excellently agreed (1.3 and 1.8% average relative 
deviations) with those of this work, we did not 
show them in the figure. We formulated a new 
correlation equation for the Henry’s constants by 
using data of Sato et al. [4] and this work for PS 
and using the data of this work for PVAc. In the 
formulation, a least-squares method was used to 
determine a linear relationship between In(1/Kp) 
and (T7./T). 

CO, + PS system. 


| | _ 6.400 2.537 ry 
i eg ‘ + Z. T 


P 


CO, + PVAc system. 


l y iY, 
5B po 8 
In 6.914 + 2.924 ( (8) 


P 
These equations could represent the experimen- 
tal Henry’s constants to within average relative 
deviations of 2.5% for PS and 1.0% for PVAc. 
Experimental solubility data were correlated 
with the Sanchez and Lacombe equation of state 
(S—L EOS) [19,20]: 
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Fig. 4. Solubility of carbon dioxide in PS. 


P= - 9° =F] ima —p)+(1-2)p | (9) 


a ee MP* 
T — T= 

| ie RT*p* 
where characteristic parameters, P*, p*, and T* 
of the S—L EOS for mixture were evaluated with 


the following mixing rules. 


(10) 


P* =) b:O;P5 (11) 
Pi, = (1 —k,)(PTPT)” (12) 


(13) 


(14) 


P*/T* 


»(b;P7/TH) 
J 
ca 
_ _/p7) (16) 
20Wj/07) 
j 
In Eqs. (11)—(16), 7*, P*, p*, and r° refer to the 


characteristic parameters of component 7 in the 
pure state and k,, is a binary interaction parameter 
determined by fitting the equation to the experi- 
mental data. The pure component parameters 
used are given in Table 3. In the calculation of the 
solubilities, it was assumed that the polymer was 
monodisperse and that it did not dissolve in the 
vapor phase. 

The binary interaction parameter, k,, in Eq. 
(12) was determined so as to minimize the relative 
deviations between experimental and calculated 
solubilities at each temperature. Correlation re- 
sults are shown as solid lines in Figs. 3 and 4. The 
Sanchez—Lacombe EOS could correlate solubili- 
ties to within an average relative deviation of 
3.6% for PVAc and 1.6% for PS. The k,, varied 
almost linearly with temperature. In both the 
systems, the deviation in the correlation became 
larger at the higher solubilities. This deviation 
may be attributed to two reasons. First, S,, cor- 
rection might be overestimated. Second, the S—L 
EOS and its mixing rules might be not suitable for 


these systems. 





~; 
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Table 2 
Experimental solubility of carbon dioxide in PS 





Temperature Pressure 


(K) (MPa) (g-gas/kg-polym) 


Solubility without S,, correction 


Solubility with S,, correction Swelling S,, (S—L 
(g-gas/kg-polym) EOS) 





373.15 2.068 12.31 
2.627 15.42 
4.093 23.50 
4.629 26.38 
6.038 33.39 
6.502 35.54 
8.054 42.04 
8.585 44.27 

12.105 54.35 
12.631 55.33 
16.059 59.56 
20.067 58.47 
2.159 8.72 
2.689 10.80 
4.073 15.93 
4.732 18.23 
6.049 22.69 
6.530 24.32 
8.071 28.98 
8.566 30.41 
12.046 38.52 
12.593 39.59 
16.143 44.29 
2.167 6.88 
2.884 9.02 
4.232 12.9] 
4.795 14.45 
6.108 17.99 
6.685 19.48 
8.027 ya 
8.813 24.61 
12.165 31.13 
12.659 32.08 
16.040 36.71 
20.151 39.88 


12.74 0.014 
16.13 0.018 
25.29 0.029 
28.71 0.032 
37.53 0.042 
40.40 0.045 
49.83 0.056 
33.29 0.060 
74.16 0.083 
77.21 0.087 
98.68 0.110 
121.84 0.136 


9.04 0.011 
11.30 0.014 


17.09 0.022 
19.81] 0.025 
25.32 0.032 
27.42 0.034 
33.80 0.042 
35.88 0.045 
49.63 0.062 
51.79 0.064 
64.46 0.080 

7.14 0.010 
9.48 0.013 
13.90 0.020 
1x73 0.022 
20.09 0.028 
21.99 0.031 
26.41 0.037 
29.05 0.041 


39.60 0.055 
41.29 0.057 


51.36 0.071 
62.31 0.085 








3.2. Diffusion coefficient 


In Fig. 6, the relationship between the relative 
amount of dissolved CO, in PVAc and elapsed 
time at 313.15 K is shown as an example. The 
solid line in the figure denotes calculation results 
of Eq. (6). For this case, a diffusion coefficient of 
1.09 x 10~'° m?/s was obtained. 

Mutual diffusion coefficients of CO, in PVAc 
at 313.15 K are shown in Fig. 7 and Table 4. The 
diffusion coefficient exhibited strong concentra- 
tion dependence and increased with increasing 


solubility. Literature data of Sasaki et al. [27] 
obtained by using a pressure decay method also 
shown in the figure. Although the diffusion coeffi- 
cients of this work agreed with the literature data 
in the lower concentration region, our measure- 
ments in the higher concentration region were an 
order of magnitude higher than those of Sasaki et 
al. [27]. 

Uncertainty of the diffusivity for CO, + PVAc 
system at 313.15 K was estimated to be less than 
28%, which includes an assumption of 5% in the 
S,, uncertainty. It was found that uncertainty of 
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Fig. 5. Henry’s constant for carbon dioxide in PS and PVAc. 


polymer thickness, L caused a large uncertainty in 
the diffusion coefficient. 

Mutual diffusion coefficients of CO, in PS are 
shown in Fig. 8 and Table 5. The mutual diffu- 
sion coefficients of CO, in PS were correlated with 
the free volume model of Kulkarni and Stern [28]. 
The self-diffusion coefficient D,.., is given by the 
expression [29]: 


Dyap = RTA; exp| aa (17) 
f 

where R is gas constant, 7 the absolute tempera- 

ture, and A, and B, are characteristic parameters. 

The v; is volume fraction of the free volume of the 

system comprised of gas and polymer was defined 

as follows [28]: 


vf T, F, P,) 


= v,,(T,, P., 0) + «(T — T,) — B(P — P.) + yd; 
(18) 


where v,;, is the free volume fraction at the glass 


Table 3 
Characteristic parameters for Sanchez~-Lacombe EOS 





CO, + PVAc 


313.15 K 


Oo ex. 
——— Eq. (6) 


[W,(f) - W,(0)] / [W,() - W(0)} (-) 











Time x10°* (s) 


Fig. 6. Example of sorption curve for carbon dioxide and 
PVAc system pressure, 3.589—3.979 MPa; solubility, 139.0— 
158.2 g-gas/kg-polym; diffusion coefficient, 1.09 x 10~ !° m?/s. 


transition temperature and atmospheric pressure, 
a the thermal expansion coefficient, 6 the com- 
pressibility y the concentration coefficient of the 
free volume, T, the glass transition temperature, 
P, the standard pressure (= 0.1 MPa), and @, is 
the volume fraction of gas. The value for v,, was 
fixed at 0.025 for amorphous polymers and « was 
fixed at 4.8x10~* per K. These values were 
investigated by Williams et al. [30] in the tempera- 
ture range of T,< T7<T,+ 100. The compress- 
ibility was determined from the following 
equation [28]: 


l 
- ef = (19) 


l—e V 

where f’ is the isothermal compressibility of the 
pure polymer obtained from PVT data [18]. The 
v* in Eq. (19) is the free volume of polymer at 


temperature, T and pressure, P, and was calcu- 
lated from the following equation. 





Substance P* (MPa) p* (kg/m?) 


T* (K) References 








Carbon dioxide 720.3 1580 
PVAc 504.2 1282 
PS 387.0 1108 


208.9 +.0.459T—7.56 x 10-47? [26] 
592.0 [17] 
739.9 [4] 








Y. Sato et al. /J. of Supercritical Fluids 19 (2001) 187-198 





T Ul T T T 
CO, + PVAc 
313.15K 





Mutual Diffusion Coefficient (m2/s) 


O. This work 
— Sasaki et al. (1991) 





A | 1 1 1 rn , a | rn 


0.1 0.2 0.3 








Solubility (g-gas/g-polym.) 
Fig. 7. Mutual diffusion coefficient of carbon dioxide in PVAc. 


vt = v,, + a(T — T,) + BP, 


vg + (T— T,) + B'P, (20) 
wi 1+ p’'P. 





The relationship between self-diffusion coeffi- 


cient D,.., and mutual-diffusion coefficient D,,., 
was obtained from [31]: 
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Fig. 8. Mutual diffusion coefficient of carbon dioxide in PS. 


where x; is mole fraction of i, 4, was chemical 
potential of gas in polymer. Mn was used for the 
calculations of x;. The chemical potential gradient 
of the gas in Eq. (21) was obtained from numeri- 
cal calculation using the SL EOS. To develop a 
correlation for the mutual diffusion coefficients 
(Eq. (21)), the characteristic constants A, and Bg, 
in Eq. (17) and y in Eq. (18) were used as fitting 
parameters to represent the experimental diffusion 
coefficients. The fitting parameters determined 
were A,=1.004x10~'? m?*mol/(sJ), By = 


Experimental diffusion coefficient of carbon dioxide in PVAc at 313.15 K 








Pressure range Average pressure 
(MPa) (MPa) 


Average solubility (g-gas/kg-polym.) 


Diffusion coefficient of CO, x 10!? (m?/s) 








0.199-0.459 0.329 10.08 
0.459-0.858 0.659 21.43 
0.858-1.271 1.065 35.48 
1.271-1.949 1.610 55.49 
1.949-2.375 2.162 76.72 
2.375-2.818 2.597 94.55 
2.818—3.260 3.039 113.7 
3.260—3.589 3.425 131.3 
3.589-3.979 3.784 148.6 
3.979-4.449 4.214 170.5 
4.449-4.978 4.714 197.7 
4.978—5.407 5.193 225.6 
5.407-5.916 5.662 255.5 
5.916-6.375 6.146 289.0 
6.375-6.876 6.626 326.3 


2.15 
3.76 
6.81 
14.7 
25.9 
39.4 
60.9 
81.3 
109 
131 
182 
223 
279 
366 
554 
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Table 5 
Experimental diffusion coefficient of carbon dioxide in PS 








Average solubility Diffusion coefficient of CO, x 10'° 
(g-gas/kg-polym) (m/s) 


Temperature Pressure range Average pressure 
(K) (MPa) (MPa) 








373.15 2.068—2.627 2.348 14.44 0.81 
4.093-4.629 4.361 27.00 1.14 
6.038—6.502 6.270 38.97 1.46 
8.054-8.585 8.320 51.56 1.67 
2.159-2.689 2.424 10.17 3.01 
4.073-4.732 4.403 18.45 4.72 
6.049—6.530 6.290 26.37 4.42 
8.071-8.566 8.319 34.84 5.33 
2.167-2.884 2.526 8.31 9.24 
4.232-4.795 4.514 14.82 10.5 
6.108—6.685 6.397 21.04 10.0 
8.027-8.813 8.420 27.73 9.90 








0.1017, y = 0.6895. Correlation results of the mu- 
tual diffusion coefficients are shown by the solid 
lines in Fig. 8. The free volume theory could 
correlate the diffusion coefficients to within 11% 
average relative deviation. 


4. Conclusions 


Solubilities of diffusion coefficients of carbon 
dioxide in PS and PVAc were measured with a 
new apparatus using a MSB technique. Experi- 
mental solubility data in PVAc and PS were cor- 
related with the Sanchez and Lacombe equation 
of state to within an average relative deviation of 
3.6 and 1.6%, respectively. Experimental diffusion 
coefficients in PS were correlated with the free 
volume model of Kulkarni and Stern to within an 
average relative deviation of 11%. 
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Abstract 





The polystyrene//CO,/ethanol system has been investigated in the supercritical domain using infrared absorption 
(FTIR) and Raman scattering spectroscopies. The pressure and temperature were fixed in the range 20-35 MPa at 
423 K for CO,/ethanol mixtures with proportion in weight from (70/30) to (58/42). These conditions are those used 
in Our previous study concerning the fractionation of this macromolecule. A surprising finding reported before was 
that the CO,/ethanol mixture made of two non-solvents of PS at ambient conditions is able to dissolve PS chains 
under high temperature and pressure. In the present article, this question will be addressed at a microscopic level from 
a discussion of the intermolecular interactions between the polymer and the medium and between the CO, and 
ethanol components provided by the analysis of the spectra. © 2001 Elsevier Science B.V. All rights reserved. 


Keywords: Polystyrene; CO,/ethanol mixtures; Raman spectroscopy 











1. Introduction 


The fractionation and the purification of poly- 
mers using supercritical fluids (SCF) is a well 
established approach in material science process- 
ing [1—4]. In this context, we have recently ap- 
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plied this technique to fractionate two chemically 
different polymers, namely poly(ethylene oxide) 
(PEO) and polystyrene (PS), with a supercritical 
mixture made of carbon dioxide (CO,) and etha- 
nol [5,6]. An interesting result came out from the 
study of PS. Indeed, both CO, and ethanol are 
non-solvents of this macromolecule at ambient 
conditions (T= 298 K; P=0.1 MPa). However, 
although at high pressure and temperature CO, 
remains a non-solvent, the addition of amount of 
ethanol allows to dissolve PS chains. More pre- 
cisely, we have shown that, at 35 MPa and 423 K, 
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a binary mixture CO,/ethanol with proportions in 
weight ranging from (70/30) to (58/42) is able to 
solvate chains having molar masses, which are, 
respectively, 3800 and 8000 g mol~'. The results 
obtained so far in this study have clearly provided 
an insight on the solvation process at the macro- 
scopic level. However, this study has also raised 
an issue concerning the nature of the interactions 
involved in the complex PS//CO,/ethanol system 
in order to have a better understanding of the 
solvation process on a microscopic scale. 

It is noteworthy that studies in this research 
area are scarce. For instance, in vibrational spec- 
troscopy, only few works based upon infrared 
absorption (FTIR) have been undertaken to study 
the solvent CO,/alcohol (generally, the percent in 
weight of alcohol does not exceed 5%) under 
supercritical conditions. The aim was to measure 
the degree of intermolecular hydrogen bonding 
between methyl alcohol-d in supercritical CO,. In 
this fluid, the equilibrium between the free 
monomer and hydrogen bonded oligomers 
(mostly tetramers) has been assessed under pres- 
sure and temperature conditions ranging from 30 
to 40 MPa and from 313 to 353 K. In addition, a 
weak complex between CO, and methanol-d has 
been put in evidence [7]. Another FTIR study has 
reinforced this point of view concerning the exis- 
tence of an interaction between CO, and 
methanol [8]. 

The goal of this article is to perform a study of 
the PS//CO,/ethanol system using Raman and 
infrared spectroscopies. The first section will be 
devoted to the study of the pure components 
(CO,, ethanol) and to the presentation of the 
results concerning the characterisation of the 
CO,/ethanol mixtures. The experimental condi- 
tions have been fixed to those used in our previ- 
ous fractionation experiments, namely at 423 K, 
in the pressure range 20-35 MPa and with (58/42) 
and (70/30) CO,/ethanol weight proportions. In 
the second section, we will present a study of PS 
in the binary mixture, under the same experimen- 
tal conditions. In a preliminary work, we have 
investigated solutions of PS in a series of solvents 
(tetrahydrofuran, cyclohexane, carbon tetrachlo- 
ride...) at ambient temperature and pressure, in 
order to identify the vibrations of PS, which are 


the more sensitive probes of the solvation. Fi- 
nally, we will discuss in the last section the solva- 
tion process on the ground of the interactions 
existing between the polymer chains and its envi- 
ronment and those between ethanol and CO.,,. 


2. Experimental details 
2.1. Materials 


CO, is provided by Air Liquide. Ethanol is a 
Prolabo product (99.9% purity) and other organic 
solvents originates from Aldrich (purity, > 99%). 
Critical temperatures and pressures of CO,/etha- 
nol mixtures have been determined by Kiran and 
co-workers [9]. 

The PS sample comes from Polymer Expert. It 
is characterised by a broad molar mass distribu- 
tion, M, = 5000 g mol~!; M,, = 8800 g mol~!; 
I= 1.7 (where M,, is the number average molar 
mass; M,, the average molar mass in weight; and 
I the molar mass distribution). 


2.1.1. Mid-infrared spectroscopy 

The measurements were performed on a Biorad 
interferometer (type FTS-60A) equipped with a 
globar source, a KBr beam splitter and a DTGS 
detector. Single beam spectra recorded in the 
spectral range 400-6000 cm~!' with a 2 cm™! 
resolution were obtained by the Fourier transfor- 
mation of 100 accumulated interferograms. 

We have used a special inconel cell, which has 
four cylindrical sapphire windows, two for the 
infrared absorption measurements with a path- 
length of 12 mm (which can be varied from 3 to 
25 mm) and the others for direct observation of 
the solution using a video camera to insure that 
there is no demixing in the sample. The optical 
cell has been shown in detail in [10]. In order to 
be sure that the solution was at the equilibrium, 
mixing was performed with a stirrer put in the cell 
and the acquisition of the spectra was performed 
10 min after a pressure change. The sealing was 
obtained using the unsupported area principle. All 
windows were positioned on the flat surface of the 
inconel plug with a 100-um gold foil placed be- 
tween the window and the plug in order to com- 
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pensate for imperfections at the two surfaces. Flat 
copper seals were used to insure sealing between 
the plug and the cell body. The heating of the cell 
was performed by four cartridge heaters, which 
were disposed in the body of the cell in which two 
thermocouples were placed, the first one was lo- 
cated close to one cartridge to achieve the temper- 
ature control and the second was kept close to the 
sample area in order to insure a good temperature 
regulation with an accuracy of about AT + 2°C. 
The cell was connected via a long stainless steel 
capillary to an hydraulic pressurising system, 
which permits to control the pressure up to 1000 
bar with an absolute uncertainty of +1 bar plus 
a R.E. of +0.3%. The lower part of the cell was 
loaded with an excess of polystyrene powder over 
that needed for saturation and the cell was evacu- 
ated under vacuum. Then, it was filled with etha- 
nol at the requested weight percent concentration. 
The cell was heated up to 423 K and CO, was 
then added up to the desired pressure. 


2.1.2. Raman spectroscopy 

The Raman spectra were collected on a DILOR 
Omars 89 double monochromator spectrograph 
equipped with a charged coupled device (CCD) 
detector. The source was a Spectra Physics argon 
ion laser operating at a wavelength of 514.5 nm 
with a power of 600 mW. The polarised /,,(7) 
spectra were recorded in the spectral range 3400— 
3800 cm~' with a resolution of 2 cm~! using the 








Scattering intensity (a.u.) 
. =. = ae Se 
eee. ee 














1300 1350 
Wavenumber (cm) 


Fig. 1. Comparison of the Raman spectra of CO, in the 
CO,/ethanol mixture ((58/42), — —) and in pure CO, (—) at 
423 K and 20 MPa. 


standard 90° scattering geometry. The Raman cell 
used here was made of titanium and had the same 
geometry and sealing principle as the infrared cell. 
However, in this case, the entrance and exit win- 
dows were made of fused silica in order to per- 
form proper spectroscopic polarisation 
observation. 


3. Results and discussion 


3.1. Investigation of the supercritical CO,/ethanol 
mixture 


3.1.1. Raman scattering 

We have performed Raman scattering measure- 
ments on the supercritical CO,/ethanol mixture at 
various ethanol concentration (0, 30 and 42 wt.%) 
at 423 K for pressures ranging between 15 and 35 
MPa. Two spectral domains are of particular 
interest here, namely, the spectral range between 
800 and 1450 cm ~', which is related to the Fermi 
resonance of the CO, internal modes (see below) 
and the 2700-3800 cm~! region associated with 
the OH stretch of the ethanol molecule. 


3.1.1.1. Spectral range 800-1450 cm~'. We have 
reported in Fig. 1 the spectra of pure CO, com- 
pared with that of the supercritical (58/42) CO,/ 
ethanol mixture at 423 K and 20 MPa. The 
Raman spectra of the pure CO, displayed two 
narrow bands observed at about 1285 and 1388 
cm~!, which were due to a Fermi resonance 
between the harmonic 2v, of the bending mode 
and the symmetric stretch v, [11]. Additional very 
weak bands were also observed, which were as- 
signed to a Fermi resonance between the modes 
2v, and v, of the isotope derivatives C'’—-O'® and 
C!?_O!® and between v, +, and 3v, [12]. These 
very weak bands will not be further considered in 
the following. The comparison of these spectra 
exhibited some noticeable differences. In the mix- 
tures, we observed that the bands were shifted to 
lower wavenumbers compared with that of pure 
CO,. Moreover, this shift increases with ethanol 
concentration ranging from 1 cm~! for the (70/ 
30) mixture up to 2.5 cm~! for the (58/42) solu- 
tion. This result might be associated with the 
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Fig. 2. Comparison of the Raman spectra of ethanol in the 
(70/30) and (58/42) CO,/ethanol mixture with pure ethanol at 
423 K and 35 MPa. For clarity, the two upper curves have 
been shifted upwards (arbitrary unit). 
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Fig. 3. Evolution of Raman spectra of ethanol as a function of 
pressure in the (58/42) CO,/ethanol mixture at 423 K. For 
clarity, the three upper curves have been shifted upwards 


(arbitrary unit). 


existence of an interaction between CO, and etha- 
nol. Finally, we found that in a CO,/ethanol 
mixture of fixed composition, when the pressure 
was increased from 15 to 35 MPa, the shifts of the 
CO, Fermi diads were less than 1 wavenumber. 
F or neat CO, when the pressure was increased in 
the same range, the shift of the Fermi diads was 
also less than 1 wavenumber. For this reason, we 
can say that there is no noticeable effect of the 
pressure on the mixture [13]. 


3200 3300 3400 3500 3600 3700 3800 


3.1.1.2. Spectral range of the OH stretching mode 
of ethanol. Fig. 2 displays the Raman spectra of 
CO,/ethanol mixtures at various ethanol concen- 
tration (0, 30 and 42 wt.%) at 35 MPa and 
150°C). The spectrum associated with the OH 
stretching mode of pure ethanol at 423 K and 35 
MPa was composed of a narrow band centred at 
about 3640 cm ~' assigned to monomeric ethanol 
together with a broad band centred at lower 
wavenumber (about 3450 cm ~ ') corresponding to 
the ethanol molecules, which were hydrogen- 
bonded aggregated [14]. This spectrum was not 
affected by a pressure variation between 10 and 
35 MPa [13]. Upon addition of CO, in the mix- 
ture, there was a strong modification of the shape 
of the spectra as it can be seen in Fig. 2 for both 
(70/30) and (58/42) mixtures. The intensity of the 
broad profile associated with alcohol oligomers 
decreased at the expenses of the narrow band 
associated with monomeric ethanol. Therefore, we 
can infer that the state of aggregation of the 
alcohol molecules is a function of the mixture 
composition and is displaced towards the forma- 
tion of monomers as the CO, concentration in- 
creases in the medium [7,15]. Moreover, we 
observe also that the narrow band associated with 
the monomeric ethanol is red shifted at decreasing 
concentration of CO, in the mixture. Indeed, the 
peak maximum was, respectively, observed at 
3637, 3647 and 3650 cm~! on going from the 
pure ethanol to the (58/42) and (70/30) CO,/etha- 
nol mixtures. Clearly, the vo, vibration of the 
‘free’ ethanol molecules probes mostly the attrac- 
tive part of the intermolecular potential due to 
both CO, and ethanol surrounding molecules. It 
is very likely that in the electrostatic interaction, 
the dipolar field of surrounding ethanol molecule 
(let’s emphasise that the ‘free’ ethanol molecule is 
not H-bonded) lead to greater perturbation of the 
vibration than the quadrupole field of the CO,. 
For a given composition (58/42), the variation 
of pressure in the 20-35 MPa range has no 
noticeable effect on the broad profile associated 
with ethanol aggregates (Fig. 3). Therefore, the 
concentration of oligomers in the mixture is not 
sensitive to pressure variations in this range. In 
contrast, the band associated with monomeric 
species was shifted towards lower frequencies 
when the pressure increased from 20 to 35 MPa 
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(Fig. 3). Attractive forces may be involved to 
explain this trend. Indeed, in a Van der Waals 
mean field approximation, it can be shown that 
the attractive forces lead to a red shift of the band 
centre, which is a linear function of the density 
(see insert in Fig. 3) [16]. Therefore, as before we 
may conclude that the monomer is in interaction 
with surrounding CO, and ethanol molecules and 
that the strength of the attractive interaction in- 
creases with the pressure [7,8,15]. 


3.1.2. Infrared absorption 

The strong absorption of the CO,/ethanol mix- 
ture precludes to obtain information about inter- 
nal modes of CO, and ethanol in the spectral 
range between 500 and 4000 cm~'!. This draw- 
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Fig. 4. Evolution of the combination bands of CO, as a 
function of pressure (20 —, 25 —, 30 — and 35 MPa —) in 
the (58/42) CO,/ethanol mixture at 423 K. 
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Fig. 5. Evolution of the combination bands of pure CO, as a 
function of pressure at 423K. 


back can be circumvented by looking at the com- 
bination modes of CO, and ethanol in the spectral 
range 4000-5000 cm~ ' although they are weaker 
than the internal modes. We have used deuterated 
ethanol (CD;CD,OD), which will allow to ob- 
serve in the PS/CO,/ethanol mixture the internal 
modes of the benzene ring of polystyrene in the 
3000-3100 cm! spectral range. We have opti- 
mised the pathlength of the cell (12 mm) to detect 
with a good accuracy the combination modes of 
CO, in the range 4700—5200 cm~!. 


3.1.2.1. The CO,/ethanol mixture. We have re- 
ported in Fig. 4 the spectra of the three combina- 
tion bands of CO, in the (58/42) CO,/ethanol 
mixture at 423 K and pressure ranging between 20 
and 35 MPa. These bands observed at about 
4840, 4960 and 5100 cm~'! are, respectively, as- 
signed to the combination modes 47, + 13, v, + 
2v,+ 3 and 2v, + v3 [11]. Clearly, it appears that 
the pressure has no appreciable effect on the 
shape and the band centre of the three bands. The 
intensity enhancement observed when the pressure 
was raised was ascribed to the increase of the 
density of the medium. These measurements have 
also been performed on pure CO, under the same 
temperature and pressure conditions. The results 
reported in Fig. 5 show that all the three profiles 
are now sensitive to the pressure variation. In- 
deed, we notice that the band shapes are modified 
and that their intensity increases when the pres- 
sure is raised. In contrast, there is no noticeable 
variation of the band centre. The variation of the 
band shape can be qualitatively explained as fol- 
lows. At low pressure, each band has the shape of 
the P—R vibrational—rotational transition envelop 
characteristic of a parallel transition of a linear 
molecule [11]. As the pressure is increased, the 
doublet structure collapses leading to a single 
profile centred at the mean frequency of the dou- 
blets. These modifications are well known in the 
spectroscopy of gaseous state [17]. Indeed, in- 
creasing the pressure, molecular collisions are 
more frequent and therefore, the rotation of CO, 
molecules are more perturbed [7,18,19]. The en- 
hancement of the intensity is related to the in- 
crease of the medium density. Finally, in order to 
estimate the influence of the ethanol addition to 
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Fig. 6. Comparison of the combination bands of CO, in the 
(58/42) CO,/ethanol mixture (— —) with that of pure CO, (—) 
at 423 K and 20 MPa. 


pure CO,, we have compared the spectrum of 
pure CO, with that of the (58/42) CO,/ethanol 
mixture under the same experimental conditions 
(Fig. 6). Clearly, the three bands observed for the 
mixture are always shifted to lower frequencies 
compared with that of pure CO,. This shift varies 
approximately from 5 cm‘ for the (70/30) CO>/ 
ethanol mixture up to 9 cm~! for the (58/42) 
solution. The main result concerns the shape of 
the profiles. Indeed, comparing these two spectra, 
we remark that the ‘doublet’ structure observed 
for pure CO, on the three bands vanishes when 
ethanol is added in solution. All these observa- 
tions allow to infer that an interaction exists 
between ethanol and CO>, which leads to an 
hindrance of the CO, molecule rotation. 


3.2. Polystyrene /solvents 


3.2.1. PS in solvents at ambient conditions 

In all the experiments performed in the follow- 
ing, we have used a polydisperse polystyrene sam- 
ple with a broad molar mass distribution (see 
Section 2.1). In order to test if the molecular 


Table | 


weight of the polystyrene sample has any influ- 
ence on the infrared and Raman spectra, we have 
performed vibrational spectroscopic measure- 
ments on a monodisperse polystyrene sample dis- 
solved in CS,. No difference was observed 
between the spectra associated with the two sam- 
ples. Such result is not puzzling. Indeed, vibra- 
tional spectroscopy is mainly sensitive to the local 
structure around a specific oscillator. 


3.2.1.1. Raman scattering. We have performed 
Raman spectroscopic measurements on PS dis- 
solved in a number of selected liquid solvents in 
view of their various solvation power, namely, in 
tetrahydrofuran (THF) (1 mol of monomer per }), 
in carbon disulfide (CS,) (10~* mol of monomer 
per 1), and in chloroform (10~7 mol of monomer 
per 1). In view of the low concentration of the 
solute, we have been able to detect only the 
symmetric breathing mode of the benzene ring of 
PS, centred at 1000 cm~ '. We observe that there 
is no dependence of the mode with the nature of 


the solvent. 


3.2.1.2. IR spectroscopy. The same methodology 
has been used here. The solvents have been cho- 
sen in order to detect a maximum number of PS 
vibrational modes. We have analysed the follow- 
ing solutions, PS//THF (10~' and 1 mol of 
monomer per 1); PS//CS, (10~7 and 10~' mol of 
monomer per 1); PS//1,4-dioxan (10~? and 107! 
mol of monomer per 1); and PS//carbon tetrachlo- 
ride (CCl,) (10~! mol of monomer per 1). Among 
the observed modes, that centred at 699.9 cm~' 
associated with the out of plane benzene ring 
deformation is slightly dependent on the nature of 
the solvent (Table 1). If we compare the band 
centre frequency with that reported for a PS film, 
we detect a shift of about 3 and 1.5 cm! to 
higher frequencies in the case of 1,4-dioxane and 


Solvent effect on the PS band (the out of plane benzene ring deformation), centred at 699.9 cm~! (PS concentration, 10~' mol of 


monomer per |) 

CS, (cm~'!) Cyclohexane (cm ~ ') CCl, (cm~') 
Re 2 ath aE eho 

698.5 698.5 699.6 

bE OER RR e0, stn OO ee Pai 


PS film (cm~') 








THE (cm~'!) 1,4-Dioxane (cm~ ') 








701.2 702.7 
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Table 2 


Solvent effect on the PS bands (CH stretch of the benzene ring) centred at 3026.4 and 3060 cm~! (PS conceniration, 10~! mol of 


monomer per 1) 





PS film (cm~') 1,4-Dioxane (cm~ ') 


CS, (cm~') 


THF (cm~') Cyclohexane (cm~') CCl, (cm~') 





3026.4 3026.4 3027.4 
3060 3060.6 3061.4 


3028.4 3028.4 3029 
3062 3063.3 3065.3 





THF, respectively, whereas a shift to lower fre- 
quencies of about 1.5 cm~' is observed in CS, or 
cyclohexane. In CCl,, the band position is the 
same than that reported for the PS film. The 
bands of PS centred at 1454, 1494 cm~! and in 
the spectral range 2900-3000 cm~ ' are hampered 
by strong absorption of the solvents. However, we 
have detected in the spectral range between 3000 
and 3100 cm~' two bands centred respectively at 
3026.4 and 3060 cm~', which are assigned to CH 
stretching modes of the benzene ring of the PS 
molecule [20]. Taking the frequency reported for 
the PS film as a reference, we observe small 
frequency shifts to higher frequencies of about 
1-3 cm~! for PS dissolved in THF, cyclohexane, 
CCL, and 1,4-dioxane and we notice that CS, does 
not influence the band position (cf. Table 2). 
Finally, the polymer concentration does not affect 
the spectroscopic observations (band centre and 
width). Using these measurements of PS dissolved 
in liquid solvent at ambient conditions, we were 
able to calculate the molar absorption coefficient 
of the band centred at 3026.4 cm~! from its 
integrated area; The value is 580 | mol~' cm~? 
and found to be independent on the nature of the 
solvent within experimental uncertainties. 


3.2.2. PS in CO,/deurated ethanol mixture 

It’s challenging to record the PS spectra in the 
CO,/ethanol mixture. Indeed, in view of the low 
solubility of the PS in the solution, the spectrum 
is very weak and obscured by the absorption of 
CO,/ethanol mixture. The only accessible spectral 
range is limited to the CH stretching of the PS 
benzene ring (3000-3100 cm~') and only if 
deuterated ethanol is used. Knowing the solubility 
of PS in the mixture as reported in previous 
studies [6], we have optimised the pathlength of 
the cell (12 mm) to detect with a good accuracy 


the CH stretching mode of the benzene ring of 
polystyrene at about 3030 cm~!. We have mea- 
sured the spectra of PS in the spectral domain 
3000 and 3100 cm ~' at 423 K and 35 MPa for PS 
in the (58/42) and (70/30) mixtures (Fig. 7). It is 
found that the band centre of the CH stretch is at 
about 3030 cm ' whatever the ethanol concentra- 
tion in the medium is. This value is in the range of 
those reported for PS dissolved in liquid solvents 
at ambient conditions (Table 2). We note also 
that it is in the CO,/ethanol mixture that the 
greater blue shift (compare with PS film) is ob- 
served. This finding is difficult to rationalise. In- 
deed, we know that the shift of the band centre 
frequencies observed in vibrational spectroscopy 
results from a subtle competition between attrac- 
tive and repulsive forces. A glance to Table 2 
shows that it’s difficult to assign the spectral shifts 
observed to a specific molecular properties. In- 
deed, it’s not obvious to correlate multipole mo- 
ment, which gives some insight about attractive 
forces as well as the molecular shape (linear, 





Perry 


Absorbance 











3010 3015 3020 3025 3030 3035 3040 


Wavenumber (cm ') 


Fig. 7. Evolution of the CH stretching vibration mode of the 
benzene ring of PS in the (70/30) and (58/42) CO,/ethanol 
mixture at 423 K and 35 MPa. The spectra of PS film is 
reported for comparison. 
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Fig. 8. Evolution of the CH stretching vibration mode of the 
benzene ring of PS in the(58/42) CO,/ethanol mixture as a 
function of pressure at 423 K. 


planar or spherical), which influences the repul- 
sive interactions. In the absence of information on 
the structure and dynamics involved in these com- 
plicated mixtures, any attempt in these analyses 
seems to us extremely questionable. We have also 
performed some measurements at constant etha- 
nol concentration as a function of pressure (see 
Fig. 8). From this set of data, it is possible to 
calculate the PS concentration in the (58/42) CO,/ 
ethanol mixture using the value of the molar 
extinction coefficient, which has been evaluated 
above. The results are reported in Table 3 and 
compared with those evaluated using UV spec- 
troscopy [6]. An overall good agreement was ob- 
served for the two sets of data (within 15% of 
errors), which ensured the validity of the present 
and previous results and showed that the current 
experimental condition to obtain a good equi- 
librium of the solutions was achieved. 


4. Conclusion 


We summarise here the main findings coming 
out from the investigation of the PS//CO,/ethanol 


Table 3 


Concentrations of PS in the (58/42) CO,/ethanol mixture at 423 K 


mixture. The study of the CO,/ethanol solvent 
reveals that ethanol is partitioned between 
monomers and H-bonded aggregates for which 
the relative concentration depends only upon the 
CO, concentration. Moreover, a weak interaction 
between the CO, and ethanol molecules has been 
put in evidence. These results corroborate experi- 
mental investigations on CO,/methanol [7,8,21] 
and CO,/ethanol [22]. The nature of the interac- 
tion might be of the Lewis acid—base type in 
which CO, carbon atom and the oxygen atom of 
the ethanol play, respectively, the role of the 
electron acceptor and donor. This conclusion 1s 
supported by ab initio theoretical work on the 
H,O-—CO, system [23,24] and from experimental 
studies on polymer impregnation with CO, [25]. 
Therefore, the PS//CO,/ethanol system is com- 
posed of several species, namely, PS chains, 
monomeric ethanol molecules interacting with 
CO,, aggregated ethanol molecules and _ finally 
CO, molecules. The ethanol is in large excess 
compared with the polymer. Estimations indicate 
that there is, respectively, about 1000 and 4000 
ethanol molecules per styrene unit in the SC 
(70/30) and (58/42) mixtures. Therefore, even if 
ethanol molecules are partly aggregated, there 1s 
still a large amount of ethanol monomers, which 
can participate to the solvation process of PS 
chains. If we turn out to PS, we have shown that 
it may locally interact with the solvent medium as 
indicated by the perturbation of the CH benzene 
vibrations. It seems that the benzene ring is likely 
to be involved in the solvation process. It is 
noteworthy from our previous works on PS frac- 
tionation that ethanol is responsible of the PS 
solvation in the supercritical CO,/ethanol mix- 
ture. Therefore, an interaction between benzene 
rings of the PS and the ethanol molecule may 
explain the solvation process. 





P (MPa) 20 


30 





C? (g monomer per g) IR 
C} (g monomer per g) UV 


1.1 x 1073 
8 x 10~* (28.8 MPa) 








P. Lalanne et al. /J. of Supercritical Fluids 19 (2001) 199-207 207 


The present investigation has also revealed the 
stringent limitations encountered in experimental 
studies of polymer diluted in supercritical binary 
mixtures using vibrational spectroscopy. Indeed, 
the low solubility of the polymer, which leads to a 
low absorption or scattering and the fact that the 
strong contribution of the solvent leaves few ac- 
cessible spectral range to the investigations of the 
solute, are very severe constraints. 

Finally, much more theoretical works remain to 
be done in order to obtain information on the 
intermolecular interaction at a microscopic level 
from the spectroscopic shift and band shapes. 
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Abstract 


BaTiO, powders from sol-gel derived gels were prepared using two different drying methods. In addition to the 
conventional drying of gels in air at 90°C the supercritical CO, drying method was also used. Results showed that 
the properties of BaTiO; powder produced by supercritical drying with CO, are better. The grain surface is less 
contaminated as a result of the supercritical drying and the microstructure development during sintering leads to a 
homogeneous fine-grained microstructure. © 2001 Elsevier Science B.V. All rights reserved. 
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1. Introduction 


An important area of investigation in the study 
of ceramic materials is related to ceramic powders 
with small grains, in particular the study of nano- 
sized powders, used during the preparation of 
electronic ceramics. Many different methods of 
synthesis for fine ceramic powders are well 
known, for example: co-precipitation, spray dry- 
ing, spray freezing, hydrothermal synthesis and 
the sol-gel procedure [1-3]. Among these pro- 
cesses, the sol-gel procedure is_ particularly 
promising since it enables the control of the 
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homogeneity and reactivity of the ceramic pow- 
ders to a large extent. Sol-gel techniques are in- 
creasingly used for the production of high 
technology ceramics, due to superior electronic 
properties, mechanical stability and chemical pu- 
rity. One of the important steps during powder 
synthesis, using the sol-gel procedure, is the dry- 
ing of the gels. Many important powder proper- 
ties, for example: the homogeneity, the size of the 
agglomerates and the extent of agglomeration are 
developed during this important step of powder 
synthesis. The drying of gels is certainly a key 
process when it comes to the final properties of 
the powder, this is particularly true when con- 
trolling the powder agglomeration. The essential 
electronic and mechanical properties, reflected in 
the microstructure, are determined to a large ex- 
tent by the drying stage of the gels. 
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PIR pressure indicator and regulator 
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Fig. 1. Experimental apparatus for supercritical CO, drying. 


Gels are usually dried for long times in air at 
elevated temperatures. In the course of our work 
we have prepared fine-grained powders by drying 
gels using supercritical CO, [4] and we have com- 
pared their properties with those of classically 
prepared powders. 


2. Experimental 


Different BaTiO, powders were prepared by 


two different drying methods for sol-gel derived 
BaTiO, gels. Besides the conventional drying of 
gels in air at 90°C (powder A) supercritical drying 
with CO, was also applied (powder B) to compare 
the properties of both BaTiO, powders. 


2.1. Sol-gel synthesis 


For the preparation of BaTiO, gels, titaniu- 
m(IV) isopropoxide, Ti[OCH(CH;),]4, (ALFA, 
purity 99.99%) and barium ethoxide, 
Ba(OCH,CH;),, (ALFA, purity 99.99%, 10 wt.% 
solution in pure ethanol) were used in equimolar 
amounts (Ba:Ti= 1:1). The synthesis of Ba—Ti 
precursor alkoxide mixtures was performed in a 
glass reactor in a dry nitrogen atmosphere and 
stirred for 48 h at room temperature. After that, 
CO,-free distilled water was injected into the pre- 
cursor solutions in excess amounts (molar ratio 
Ba:H,O = 1:25). The hydrolized solutions were 
aged for 30 days at room temperature. 


2.2. Supercritical CO, drying 


After the gelation was terminated, the gels were 


Fig. 2. TEM images of the BaTiO, powders: (a) prepared by drying of gels in air; (b) prepared by supercritical CO, drying. 
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Fig. 3. SEM images of BaTiO; powders: (a) powder A pre- 
pared by drying in air; (b) powder B prepared by supercritical 
CO, drying. 


dried in air at 90°C (powder A) and in super- 
critical CO, (powder B). 

Aerogels, with their special properties, can 
only be produced when the solvent is removed 
in the absence of capillary forces in the gel 
pores, because even small capillary forces col- 
lapse the gel structure. 

In supercritical carbon dioxide drying, alcohol 
is replaced by carbon dioxide at conditions 
above the critical conditions of CO, (Tc= 
31.1°C, Po = 7.36 MPa). All supercritical carbon 
dioxide drying experiments were performed at 
the temperature of 40°C and the pressure of 100 
bar. 


The extractor (Fig. 1) is completely filled with 
solvent (methanol) in order to minimize the 
evaporation of the solvent from the alcogel and 
to avoid cracks during pressure build up. CO, is 
pumped into the extractor to a pressure above 
its critical pressure. Before entering the extrac- 
tion, CO, is heated to a temperature above its 
critical temperature. During drying, the flow of 
CO, to the extractor is constant, 0.54 kg/h, and 
controlled independently of pressure. In the sep- 
arator, the solvent is separated from carbon 
dioxide. After 5 h, when the solvent is com- 
pletely replaced, the pressure is slowly reduced 
(AP/At = 3 bar/min) to ambient pressure. 


2.3. Characterization 


The dried gels were then characterized. The 
powders were inspected with scanning electron 
microscopy (SEM) and transmission electron mi- 
croscopy (TEM). The specific surface was mea- 
sured using the BET method and X-ray powder 
diffraction patterns were taken. The crystallite 
size of BaTiO, was estimated using the Scherrer 
equation, from the (200) diffraction peak. The 
surface of powders eroded with argon and the 
impurities concentration depth profiles were ana- 
lyzed using Auger electron spectroscopy (AES). 
The erosion rate was 1 nm/min. In addition, the 
shrinkage of the samples, prepared from both 
powders, was determined and the microstructure 
of polished and etched samples was investigated. 


3. Results and discussion 


Fig. 2 shows typical TEM images of the dried 
BaTiO; powders A and B. It can be seen that 
the synthesized powder is composed of crystal- 
lites whose size can be measured in nanometers. 

In Fig. 3 the images of dried powders A 
(dried in air) and B (dried with SC CO,) taken 
with a SEM are shown. In a comparison of 
both images one can see that powder A, which 
was dried in air, is more agglomerated. 

The X-ray powder diffraction patterns of 
dried powders and of calcined powders A and B 
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Fig. 4. X-ray diffraction patterns of synthesized BaTiO, powders; (a) powder A, dried in air, (b) calcined powder A, (c) powder B, 


dried with CO, and (d) calcined powder B. 


are shown on Fig. 4. The X-ray diffraction pat- 
tern of powder A (dried in air and calcined) shows 
only the diffraction peaks of BaTiO;, whereas the 
diffraction pattern of powder B (dried in SC CO,) 
shows also the diffraction peaks of BaCO,, which 
disappear after the calcination. 

The BaTiO; crystallite size was determined by 
the Scherrer relation: 


d,. = 0.941/B cos 0 (1) 


where 2 (1.5406 A) is the characteristic wave- 
length, f is the true half peak width and @ half the 
diffraction angle of the centroid of the peak (in 
degrees). 

The size of the crystallites obtained from the 
Scherrer equation is larger than that estimated 
from the specific surface area (S) for both pow- 
ders. The differences in the average size between 
d, and that obtained from the specific surface d,, 
(d, = 6V,,,/S, where V,,, is mole volume of BaTiO; 
and (S) is the measured specific surface) appear 
when the shapes of the crystallites can not be 
taken as perfect spheres. 

After the calcination of powders A and B at 
1100°C their specific surface decreases drastically 
from more than a hundred to some few square 
meters per gram (Tables 1 and 2). During the 


calcination the dense agglomerates are formed. In 
Fig. 5 the SEM images of calcined powders A and 
B are shown and it can be visually estimated that 
the average size of the agglomerates is around 150 
nm for both powders. 

However, when we estimate the size of the 
particles from the specific surface of powders A 


Table 1 

The specific surface area (S) and the average grain size d,. and 
d, obtained from the diffraction data and from the specific 
surface measurements, respectively 





S (m’/g) d. (nm) d,(nm) 


Powder A 136.6 15.4 7.3 
Powder B 147.5 16.6 6.7 


Sample code 








Table 2 

The specific surface area (S) of powders A and B calcined at 
1100°C and the average grain size of particles estimated from 
the Scherrer relation and from the specific surface area, respec- 
tively 





S (m?/g) d,. (nm) d, (nm) 
Powder A 6.4 50.3 155 


Powder B 4.3 17.3 160 


Sample code 
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Fig. 5. SEM images of powders calcined at 1100°C for 2 h: (a) 
powder A obtained by drying of gels in air; (b) powder B dried 
by supercritical CO . 


and B (Table 2), a similar value to the average 
grain size determined from the SEM images of 
calcined powders (Fig. 5a and b) is observed. 
On the other hand, the average crystallite size 
obtained from the powder diffraction data are 
17.3 nm for powder A and 50.3 nm for powder 
B. On the basis of these findings it can be con- 
cluded that the calcined powders A and B are 
composed of dense agglomerates, which are 
composed of much smaller BaTiO, crystallites. 
Both powders have a low tap density, approx- 
imately 0.97 g/cm’ for powder A and 0.45 g/cm? 
for powder B. From powders A and B, pellets 
with a raw density of 45 and 50% respectively, 


were prepared by pressing the powder with a 
pressure of about 10° Pa. Both powders exhib- 
ited bad compressibility due to a low tap den- 
sity. 

In Fig. 6 the shrinkage and the shrinkage rate 
spectra of samples A and B during sintering are 
shown. Sample A shows the maximum shrinkage 
rate at 1200°C and a final shrinkage of 40%, 
while sample B shows a maximum shrinkage 
rate at 1229°C and a final shrinkage of 37%. 

The final density of the samples after sintering 
in air, with respect to the theoretical density, 
was 95% for sample A and 91% for sample B. 
The microstructures of the samples after sinter- 
ing are shown in Fig. 7a and b. The samples 
show a relatively dense microstructure with an 
average grain size of 40-50 wm for samples A 
and 10—20 wm for samples B. The difference in 
the final grain size is associated with the initial 
powder properties and could be explained in 
terms of the powder preparation procedure. 

For BaTiO, discontinuous grain growth is the 
main process of microstructure development. In 
this process, a small fraction of the grains in a 
matrix of similar grains outgrow on account of 
neighboring grains. For this discontinuous grain 
growth to occur, firstly nucleus grains must be 
formed, which then grow much faster in com- 
parison to the other grains in the matrix [5]. The 
final microstructure depends to a large extent on 
the number of these nucleus grains. 

In the case of powder A, the impurities, i.e. 
the disintegration products formed during gela- 
tion, are deposited on the grain surface of pow- 
ders A before the sintering. Because of this, the 
grain growth during sintering is retarded since 
the mass transport between the neighboring 
grains is hindered. This has a consequence that 
less grains outgrow in the larger nucleus grains. 
The smaller number of large grains, which might 
act as nucleus grains and grow on account of 
other neighboring grains, has the consequence 
that the final grain size is large, Fig. 7a. In the 
case of powder B, the grain surface is less con- 
taminated as a result of drying in supercritical 
carbon dioxide. CO, dissolves and washes away 
the organic disintegration products formed dur- 
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ing gelation, which would otherwise retard the 
grain growth and/or the formation of the nucleus 
grains. A greater number of larger grains that 
might outgrow in the nucleus grains has the con- 
sequence that the final grain size is smaller, Fig. 
Tb. 

The investigation of the grain surfaces of pow- 
ders A and B by Auger electron microscopy 


(AEM) confirms such an explanation. Powder A 
contains more carbon rich impurities on the grain 
surfaces, which are the remains of the disintegra- 
tion products formed during the gelation and 
drying, Fig. 8a. The amount of carbon rich disin- 
tegration products on the grain surfaces of pow- 
der B is smaller Fig. 8b, i.e. the concentration 
profile of carbon vs. the ionic erosion is smaller. 
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Fig. 6. The shrinkage spectra of powders: (a) powder A obtained by drying of gels in air; (b) powder B dried by supercritical CO). 
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Fig. 7. Microstructure of BaTiO, powders: (a) dried in air; (b) 
dried by supercritical CO). 


4. Conclusions 


The investigations have shown that there is a 
link between the method of drying the gels, the 
powder properties and the microstructure devel- 
opment. The differences in the morphology of the 
powders due to the preparation technique influ- 
ence the microstructure development during the 
sintering of BaTiO; powders. The application of 
supercritical CO, during the treatment of the gels 
of BaTiO; was found to be effective. However, 
CO, reacts with unreacted BaO and forms ingre- 
dients of BaCO,. Thus, the use of some other 
supercritical fluid, which does not react with BaO, 
would in this case be more convenient. 
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